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ABSTRACT 

Post-starburst, or E+A galaxies, are the best candidates for galaxies in transition from being gas- 
rich and star-forming to gas-poor and passively-evolving as a result of galaxy-galaxy interactions. To 
focus on what E-t-A galaxies become after their young stellar populations fade away, we present the 
detailed morphologies of 21 E-l-A galaxies using high resolution HST/ACS and WFPC2 images. Most 
of these galaxies lie in the field, well outside of rich clusters, and at least 11 (55%) have dramatic tidal 
features indicative of mergers. Our sample includes one binary E-l-A system, in which both E-|-As 
are tidally disturbed and interacting with each other. Our E-|-As are similar to early types in that 
they have large bulge-to-total hght ratios (median B/T = 0.59), high Sersic indices, (n > 4), and 
high concentration indices (C > 4.3), but they have considerably larger asymmetry indices {A > 0.04) 
than ellipticals, presumably due to the disturbances within a few caused by the starburst and/or 
the galaxy-galaxy interaction. We conclude that E-|-As will be morphologically classified as early- 
type galaxies once these disturbances and the low surface brightness tidal features fade. The color 
morphologies are diverse, including six E-|-As with compact (0.4-1.4 kpc) blue cores, which might 
be local analogs of high-z ellipticals with blue-cores. The large fraction (70%) of E-t-As with positive 
color gradients indicates that the young stellar populations are more concentrated than the old. These 
positive color gradients (i.e., bluer nuclei) could evolve into the negative gradients typical in E/SOs 
if the central parts of these galaxies are metal enhanced. Our E-|-As stand apart from the E/SOs in 
the edge-on projection of the Fundamental Plane (FP), implying that their stellar populations differ 
from those of E/SOs and that E-|-As have, on average, a M/L that is 3.8 times smaller. The tilt of the 
E-l-A FP indicates that the variation among their stellar populations is closely tied to the structural 
parameters, i.e., E-j-As follow their own scaling relationships such that smaller or less massive galaxies 
have smaller M/L. We find a population of unresolved compact sources in nine E-|-As (45%), all of 
which have merger signatures. In the four E-|-As with suitable color data, the compact sources have 
colors and luminosities consistent with newly-formed star clusters. The bright end of the cluster LF 
is fainter in redder E-|-A's, suggesting that the young star clusters fade or are disrupted as the merger 
remnant ages. In summary, the morphologies, color profiles, scaling relations, and cluster populations 
are all consistent with E+As evolving ultimately into early-types, making the study of E-f-As critical 
to understanding the origin of the red sequence of galaxies. 

Subject headings: galaxies: evolution — galaxies: interactions — galaxies: starburst — galaxies: star 
clusters — galaxies: stellar content 



1. INTRODUCTION 

If some galaxies evolve from star-forming, gas-rich, 
disk-dominated galaxies (late- types) into quiescent, gas- 
poor, spheroid-dominated galaxies (early- types), we 
should find objects caught in the midst of this transfor- 
mation. The best candidates are the so-called "E-l-A", 
"K-l-A", or "post-starburst" galaxies (jPressler &: GunnI 
Il983t ICouch fc Sharpleill987l ) due to their combination 
of late- and early-type characteristics, including both a 
significant young stellar population (age < 1 Gyr) and 
a lack of on-going star formation. These galaxies have 
been spectroscopically identified by their strong Balmer 
absorption lines and absence of emission lines (e.g., [O II] 
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and Hq ) i n various environments and at a ll redshifts 
(Zabludoff ct al. 1996; Poggianti ct al. 1999; Goto et all 
12003; Blake et al.,,2004: ,Tran et al.,,2003. .20041) . 

While the cause of the abrupt end of their star for- 
mation is poorly understood, there is strong evidence 
that galaxy-galaxy tidal interactions or mergers trigger 
the starburst in many cases. First, most E-l-A galax- 
ies reside in low-density environments, such as poor 
roups, that are similar to those of star-forming galaxies 
Zabludoff et al.''1996'; 'Ouintero et al."2004'; 'Blake et al' 
^04; Balogh et al. 2005; Goto 2005; Hogg ct al. 200 ' 
Yan et al.l l2008l ). Therefore, many E-f As must arise 
from a process common in the field, such as galaxy- 
galaxy interactions, instead of a mechanism limited to 
dense r, hotter e nvironments, such as ram pressure strip- 
ping (jGunn &: G ott 1972) or strangulation (Balog h et al.l 
|2000[ ) . Secon d, a significant fracti o n of E-|-As have 



tidal features (IZabludoff et al.l [1991 lYang et al.l 12004 
iBlake et al.l 120041: iTran et all 120031 120041: iGotd l2005f) . 



2 



Yang et al. 



Third, optical and NIR colors show that the spectral sig- 
natures of E+As require enhanced recent star formation, 
rather than simply a t runcation of star f ormation in a 
normal spiral galaxies (jBalogh et al.l [20051 ). 

What will E+A galaxies become? In general, E+ As are 
bulge - dominated, high ly-concentrate d (£)uintcro et al.l 
2001: iTran et al.l [Ml iBlake et all [200l iGotol 12005 



Balogh et al. 2005), re latively gas-poor (jChang et al.l 
.200 1: Bu vle et al.]2006l ). and kinematically hot systems 
([Norton et al.ll2001h . Therefore, in a statistical sense, 
E-|-As are likely to become E/SO galaxies. However, 
due to a lack of spatial resolution in previous studies, 
we do not know whether the detailed properties of indi- 
vidual E-t-As, e.g., their bulge fractions, color gradients, 
internal kinematics, and newly formed stellar clusters, 
are consistent with their presumed e volution int o early 
type galaxies. Using HST/WFPC2, lYang et all ((200l ) 
showed that their morphological features are consistent 
with a transition from late to early types, but their sam- 
ple contained only the five bluest E+As galaxies fr om the 
Las C ampanas Redshift Survey (LCRS; Zabludo ff et al.l 
(|1996l )) and thus was not representative. As a result, we 
still do not know whether the entire population of E+As 
will evolve into E/SOs, whether there is a distinguishable 
subclass of E+As that evolves into E/SGs, or whether 
E+As evolve into typical E/SGs. 

To answer these questions, we must understand how 
well E+As match the full range of E/SO properties. Most 
fundamentally, are the global morphologies (e.g., bulge- 
to-total light ratios and concentration) of the whole 
LCRS E+A sample consistent with those of E/SOs? 
Second, E/SGs in the local universe become redder to- 
ward their center; these negati ve color gradients o rigi- 
nate from metallicity gradients (jPeletier et al.|[T990l ). In 
contr ast, E+As exhibit a wide range of color morpholo- 
gies (jYang et al . 2GG4; Yamauchi ct al. 2GG5). Can the 
color profiles of E+As evolve into those of the typical 
early types? Third, the number of globular clusters per 
unit luminosity is higher in earl y types than in late types 
([Harris fc van den Berghlll981f ). If E+As are in transi- 
tion from late to early types, one should find new star 
clusters formed during the star burst. Are there such clus- 
ters, and, if so, do their colors and numbers coincide with 
the expected evolution of the globular cluster systems of 
present-day E/SOs? Fourth, early- type galaxies lie on the 
Fundamental Plane (FP), an empirical scaling relation 
between the effective radius, the central velocity disper- 
sion, and the mean surface brightness, with remarkabl y 
small scatter ([Diorgovski fc Davis! [iQSTi : lDressleijll987D . 
Will E+As lie on the same FP once they evolve? 

To determine whether E+As evolve into objects that 
are indistinguishable from the bulk of E/SOs, we present 
HST/ACS o bservations of the 15 rem aining E+A galax- 
ies from the iZabludoff et al.l (|1996D sample. We com- 
bine these with t he previous HST/ WFPC2 observations 
of five blue E+As (lYang et al.ll2004| ) and o f one serendipi- 
tously discovered E+A (Yang ct al. 2GG6). The resulting 
high resolution imaging of the 21 confirmed E+A galaxies 
in the LCRS sample enables us to study the detailed color 
morphologies and the properties of the newly formed star 
cluster candidates at the sub-kpc scale. The detailed 
morphologies, color profiles, and cluster populations of 
many E+As are consistent with the galaxy-galaxy inter- 
action scenario. We discuss how these features will evolve 



and then relate this evolution to the properties of E /SOs. 
Using existing kinematic data ([Norton et al.l[200lh and 
our HST photometry, we also address whether or not 
the various scaling relations of E+As are consistent with 
those of E/SGs. 

This paper is organized as follows. We describe our 
E+A galaxy sample and the HST/ACS data reduction 
in fj2| In fj3l we examine the general morphology of 
E+As, including a discussion of tidal features (i 33.ip . 
surface brightness profiles (^ !j3.2p . structural parameters 
( ^3.3|) . and concentration/asymmetry measures fi i3.4p . 
The color profile s, including a clas s of E+As with lumi- 
nous blue cores ([Yang et al.ll2GG6l ). are presented in 21 
We compare the scaling relations of E+As with those of 
E/SOs in Sjni We present the properties of the newly- 
formed young star clusters in f|6l We summarize in ^J7l 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. LCRS E+A Sample 

Our HST imaging sample consists of 21 nearby E+A 
galaxies spectroscopically identified from the 11,113 
galaxy s pectra of the Las Ca mpanas Redshift Survey 
(LCRS; IShectman et al.l [loo s') with redshifts between 
0.07 and 0.18 ([Zabludoff et al. 1996) . These E+As are 
required to have strong Balmer absorption features (av- 
erage equivalent width (HPjS) > 5. 5 A) and little if any 
[O II] emission (EW[0 II] < 2.5 A). Three-quarters of 
this sample are in the field, well outside rich cluster en- 
vironments. The number of e ach E +A (e.g., EAOl) is 
adopted from IZabludoff et"an ([1996D and increases with 
increasing 4GGG A break (I?4ooo) strength. 154000 is re- 
lated to the galaxy's color — bluer galaxies have smaller 
£'40oo- For the remainder of this paper, we refer to each 
galaxy by its assigned number. 

There have been two changes to the original sample 
since published. One of the original galaxies (EA2G) was 
misclassified as an E+A due to noise in the reg ion of 
the spectral line diagnostics ([Norton et al.[ [200 1[ ) . On 
the other hand, the morphologically disturbed compan- 
ion galaxy of EAOl fies at the same re dshift and also has 
an E+A spectrum ([Yang et al.[[2GG6l ). We include it in 
our sample and refer to it as EAOIB, and to the orig- 
inal EAOl as EAGIA. The EAOIAB system is the first 
known binary E+A system and provides additional ev- 
idence that the E+A phase of galaxy evolution can be 
triggered by galaxy-galaxy interactions. 

This LCRS sample is the most extensively studied 
E+A sample, with d ata that includes rn easurements of 
internal kinematics ([Norton et aO [2001[ ). H I content 
([Chang et al.l[2001[: Wnyle et al.l[2GG6fl. radio-cont inuum 
emission (IChang et al."2001'; 'Miller &_Owen|[200T), NIR 
fiuxes and morpho logies ( Galaz 2000) , HST optical mor- 
phologies ([Yang et a l. 2GG4), and spectral diagnostics of 
nuclear activitv ([Yang et al.ll2006D . Table [H summarizes 
the basic data and inferred properties from the vari- 
ous studies. Throughout this paper, we adopt Hq = 
70 km s-i Mpc-\ Qm = G.3, and = 0.7. 

2.2. HST ACS Observations 

By combining the observations of six E+As with 
HST/ WFPC2, EAGIAB through EAG5 ([Yang et all 
[200l [2006,), with the new ACS/WFC imaging of the 
remaining 15 E+As, we now have data for the complete 
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Fig. 1. — (Left) High-contrast R band (i?702) images show the low surface brightness tidal features. (Middle) R band images for the 
WFPC2 sample (EAOIAB - EA05). (Right) Residual R baud images subtracted from the smooth symmetric model components. We 
bound each image with 4 arcsec tickmarks and include a 4 kpc horizontal scalebar. Note the diverse morphologies of E+A galaxies; tidal 
and disturbed features, dusty galaxies, blue-cores, bars, and even compact star clusters. Because EAOIA is too disturbed to be modeled 
by axisymmetric models, we restrict our analysis to EAOIB and show the residual image only for EAOIB in the top panel. 
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Fig. 2. — (Left) Same as for Fig[T] except these are ACS images of EA06-20. (Middle) Two-color composite images from the B and R 
bands. (Right) Same as for Fig[T] except for the ACS images. 




Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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ICRS sample. The ACS/WFC imaging is in the F435W 
and F625W bands, which closely match the Johnson B 
and the SDSS r bands, respectively. We chose these 
filters to match the two WFPC2 filters (F439W and 
F702W) used bv lYang et all (|2004l ). Hereafter, we wiU 
refer to magnitudes in each of these filters as i?435, -R625, 
^439, and i?702, respectively. 

The F625W band images consist of three cosmic ray 
split exposures, one short (200s) and two long (350s) ex- 
posures, while the F435W band images consist of only 
two exposures. The exposure time for each split in 
F435W ranges from 500 to 550 sec depending on the 
galaxy's location on the sky. These uneven exposure 
times per pointing and the number of splits were adopted 
to obtain both -B435 and i?625 band exposures within an 
orbit, thereby maximizing the efficiency of our HST pro- 
gram by avoiding unnecessary readout overheads. 

The images are bias-subtracted and flat-fielded with 
the standard ACS calibration pipeline. We determine 
the offsets and rotations between the exposures by cross- 
correlating the bright field stars and then "drizzle" the 
individual frames into the geometrically corrected out- 
put frame using the Pyraf/Multidrizzle package. We use 
an output pixel scale of 0.05"per pixel and the square 
interpolation kernel for the drizzle. 

We reject cosmic-rays using the standard Multidrizzle 
routine for the i?625 band images, which h ave three ex- 
posur es, and the LACOSMIC algorithms ()van DokkumI 
l2001t ) and Multidrizzle together for the ^435 band im- 
ages. Because there are only two images and each expo- 
sure is moderately deep, our i?435 band images do have 
coincident cosmic ray hits in the exposures. To address 
this problem, we identify cosmic rays in the individual 
images using LACOSMIC, fix the flagged pixels by in- 
terpolating across neighboring pixels, and then apply the 
standard drizzle process. We apply LACOSMIC in a con- 
servative way to retain real objects (e.g., star clusters) 
and visually check the rejected cosmic rays within a 20" 
radius from each E-l-A. 

We use the Vega photometric system calibrated with 
the ACS WFC zero poiii t s, 25 J79 (S435) and 25.731 
(^625) from lSirianni et all (|2005D . and cor rect for Galac- 
tic ex tinction using the reddening laws (jSirianni et al.l 
1200,51) Af625W = 2.633 E{B - V) and Af435W = 
i.l03 E{B — V) for an Sc type spectral energy distribu- 
tion (SED). The col or excess, E(B — V), at each position 
is obtained from the lSchlegel et all (|1998f l maps. Because 
the LCRS sample is at high galactic latitude, the extinc- 
tion correction is almost negligible for our sample. We 
do not include the errors from the extinction correction 
in our error budget. 

2.3. Magnitude Transformation and K-correction 

To compare our photometric measurements with our 
previous WFPC2 imaging, and with the color gradients 
(Johnson B and Cousins R) and the Fundamental Plane 
data (Johnson B and Gunn r) of early type galaxies in 
the literature, we establish the magnitude transformation 
relation between the ground, WFPC2, and ACS/WFC 
photometric systems using synthe tic photometry ofK+A 
type templates as outlined by ISirianni et al.l (|2005( ). 
Because the spectra of post-starburst galaxies depend 
strongly on burst strength and time elapsed since the 
burst, we build a library of K+A type spectra using 



stella r population synthesis models ijBruzual fc Chariot 
120031 BC03). The model stellar populations are com- 
posed of a young burst population (A-type) and an old 
underlying population (K-type) . We assume that the old 
stellar population started forming 5 or 10 Gyr ago and 
has since experienced an exponentially declining star for- 
mation rate (r = 2 Gyr). Superposed is a young stellar 
population with burst mass fraction, /burst, which repre- 
sents the fraction by mass of all stars formed in a brief, 
Atburst — 20-300 Myr, burst. We allow the burst mass 
fractions to range from 0.02 to 1.0. We also vary the 
metallicities of the old and young stellar populations be- 
tween Z = 0.004 and Z = 0.02. The purpose of this 
analysis is not to reconstruct the star formation history 
(SFH) of our galaxies, but to generate a grid of SEDs 
that covers the wide range of observed spectral parame- 
ters and broad-band colors. 

Using this grid of SEDs and the IRAF/Synphot rou- 
tine, we measure the Johnson B, ^435, ^439, Cousins 
R, Gunn r, -R702 and i?625 magnitudes over time until 2 
Gyr after the end of burst (well beyond when the galaxy 
would be classified as an E-l-A). We derive transforma- 
tion equations by fitting to the data across the model 
grid using the following linear relations, 

TOOL ^co + ci SCOL + ca SCOL^ (1) 

where TCOL and SCOL represent the target and source 
color (e.g., TCOL = S435 - i?625 and SCOL = B439 - 
iZyoa). In Table [2l we list the conversion parameters for 
the WFPC2, ACS, Landolt (i.e., Johnson B + Cousins 
R) and Thuan-Gunn systems. 

To compare the properties of E-|-As across a range of 
redshift (z = — 0.15), one should apply K-corrections. 
In principle, the X-correction can be derived from the 
known SED of the object for an arbitrary redshift and 
filter. If SEDs are not available, but there are enough 
measurements of the broadband magnitudes, one can re- 
constru ct the galaxy SED an d calculate the K-correction 
(e.g., iBlanton et al.l I2003D . However, neither flux- 
calibrated spectra nor multiple broadband colors are 
available for our sample. 

Therefore, we adopt an alternative approach using the 
synthetic photometry of the described model K-t-A SED 
library. The procedure is similar to what we use to derive 
the magnitude transformation equations. Using the grid 
of SEDs, we redshift the rest frame spectra to the appro- 
priate E+A redshift, z, measure the observed B^^^ and 
Rq25 magnitudes, and calculate the differences from the 
corresponding rest frame magnitudes (i.e., K-corrections) 
as a function of BI35— i?|25 color. For each E-l-A, we de- 
termine the K-correction at the measured (-B435— iJ^js) 
color and use the range of model values as an error esti- 
mate for the K-correction. 

3. MORPHOLOGY 

3.1. Qualitative Description 

We show images of the six WFPC2 E+A galaxies in 
Figure [TJ The leftmost panels are high contrast images 
in R and reveal low surface brightness features. The 
middle panels are the R band images (the B band images 
are far more shallow). Figure [2] shows the ACS' data 
for the remaining 15 E+As. The leftmost panels are 
the same as Figure [T] The middle panels here are the 
two-color composite images from the B and R bands. 
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Fig. 3. — E-(-A surface brightness profiles. The data are extracted using the ELLIPSE task (dark sohd line) and the model profiles 
with GALFIT (gray solid line). The dashed and dot-dashed lines are the bulge and disk components from GALFIT, respectively. For the 
barred galaxies (EA14 and 2l), we show the three components fits (Sersic + Sersic + disk; Appendix [XJ . The differences between the data 
and the model are shown in the bottom panels. The large excess of light at large radii (r ^ lOkpc) in some E+As (EA08, 10, 11, 12, 18) 
is mainly due to low surface brightness tidal features. 



We adopt a color scheme employing the arcsinh stretch 
from Lupton ct al. (2004!). We bound each image with 4 
arcsec tickmarks and include a 4 kpc horizontal scalebar. 
The right panels of both figures are the R band residual 
images resulting from the subtraction of best-fit models 
(see 



High resolution HST images enable us to identify a 
wealth of small and large scale features. Readers are 
referred to Appendix |X] for detailed descriptions of the 
morphological features. Here we briefly describe our find- 
ings. First, the morphologies of E-|-As are extremely di- 
verse, including train- wrecks, barred galaxies, blue-cores, 
and relaxed-looking disky galaxies. Given that our E-l-A 
sample was selected using uniform spectroscopic crite- 



ria, it is striking that the morphologies are so varied. 
Second, more than half (55%) of the E-|-As have tidal 
and/or disturbed features brighter than our detection 
hmit [fiji < 25.1 ± 0.5 mag arcsec"^), supporting the 
picture that galaxy-galaxy interactions/mergers are re- 
sponsible for E-l-A phase. Five E-|-As also have apparent 
companion galaxies within ^ 30 kpc that appear to be 
interacting with the E-l-A. Third, six E-|-As (30%) ex- 
hibit distinct compact blue cores with a characteristic 
size of ~ 0.5". Fourth, seven E-|-As show dust features, 
such as lanes and filamentary structures, in the two-color 
composites or the residual images. Only three of those 
have color profiles that are seriously affected by these 
dust lanes (fg]). 
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Fig. 3. — Continued. 

3.2. Surface Photometry and Model Fitting 

To compare E+A morphologies with those of early- 
and late-type galaxies in a quantitative way, we fit mod- 
els to the surface photometry. First, to obtain the sur- 
face brightness profiles (and the color profiles in we 
extract the one-dimensional azimuthally-averaged light 
profiles of the galaxies using the IRAF/ELLIPSE task. 
Because a large fraction of E-|-As are disturbed and have 
asymmetric features, we allow the model to follow the 
light distribution as closely as possible. Thus, we allow 
the center, major axis position angle, and ellipticity of 
each ellipse to change freely, but we force ellipses to not 
overlap by fitting over a limited radial range along the 
major axis. The only difference from our analysis of the 
previous WFPC2 imaging is that we forego point-spread 
function deconvolution because the errors in the drizzled 
ACS images are correlated. 

To obtain global photometric parameters, such as effec- 
tive radius r^, effective surface brightness fj,e, and Sersic 
index n, we use the two-diinensiq nal image fitting algo- 
rithm GALFIT (|Peng et al.ll2002[ ). GALFIT assumes a 
two-dimensional model profile for the galaxy with the fol- 
lowing free parameters: the {x, y) position of the center, 
Mtot (the total magnitude of the component), r^ (the 
effective radius), n (the Sersic index), q (the axis ra- 
tio defined as 6/a), the major axis position angle, and 
c (the diskiness/boxiness index, where c > indicates 
boxy). This index c plays the same role as the cos 46* 
Fourier coefficient t erm used often in isophote analysis 
(jJedrzeiewskil Il987( ) . As GALFIT explores parameter 



space, it convolves the model image with a point-spread 
function (PSF) and compares it to the data for each pa- 
rameter set. 

We use three types of model profiles with a pre- 
determined fixed sky background: 1) r^/'*-law bulge plus 
an exponential disk, 2) single de Vaucouleurs r^/*-law 
profile, and 3) single Sersic r-'^/"-law profile. Therefore, 
for each galaxy, we have three sets of structural param- 
eters, unless bulge/disk decomposition proves to be im- 
possible (e.g., EAOIA and 19). Bulge-disk decomposi- 
tion (model 1) is useful for determining whether there is 
a disk and for measuring a bulge fraction (B/T). Model 
2 is used mostly to compare the morphologies of E-|-As 
with those of normal ellipticals or SO galaxies, and we 
use the derived parameters to place E-|-As on the Funda- 
mental Plane {^S^- Model 3 is flexible enough to distin- 
guish early type galaxies (n > 3) from disk-dominated 
late type galaxies (n ^ 1). We determine the adopted 
sky value using median values from 40x40 pixel regions 
that appear free from tidal features or other large scale 
sources. The GALFIT input PSFs are evaluated using 
unsaturated stars in the ACS field with S/N higher than 
100 as identified by the IRAF/DAOPHOT package. Be- 
cause the number of available stars is often too small 
{N < 5) in each E-l-A field, we generate the time- and 
spatially-averaged PSFs using all of the stars {N ^ 90) in 
our HST observations. Wc fit the galaxy images with all 
three profiles and determine which profiles are the best 
fits based on the values and visual inspection of the 
residual images. 
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Fig. 4. — (a) Distribution of R band bulge fractions {B/T) of 15 E+As where the bulge-disk decomposition is reliable (shaded histogram). 
Except for two disky E+As (EA09 and EA17), most E+As (13/15) have B/T's greater than the median for Sa galaxy bulges (0.45). The 
median B/T of E+As (0.59) is roughly consistent with that of SO galaxies (0.63). (fe) Distribution of Sersic indices n in the R band 
for 20 E+As (except EAOIA) using a single Sersic fit (solid) and Sersic + disk fit (shaded). Most E+As (17/20) have very high Sersic 
index (n ^5), even higher than the n = 4 of de Vaucouleurs profiles of normal ellipticals, indicating that the galaxy luminosity is highly 
concentrated, (c) Sersic indices before and after masking the galaxy center. Most galaxies (65%) have a smaller n after masking, suggesting 
that the central pixels play a role in the large derived n values. 



We plot the light profiles derived from the ELLIPSE fit 
and the best fitting model profiles from GALFIT in Fig- 
ure [31 As expected, the 1-D and 2-D fits do not always 
agree, especially when there are a significant changes in 
PA or ellipticity, but generally they are consistent to 
within ^0.2 mag. Because systematic errors, such as 
that introduced by the fixed sky value, dominate, we run 
GALFIT with various plausible fitting input parameters 
and take the full range of fitting results as representative 
of the final errors. Therefore, the errors given in Table S] 
and[5]are conservative estimates of the true uncertainties. 

3.3. Light Profiles 
3.3.1. Bulge-disk Decomposition 

In Figure [4]2, we show the distribution of the R band 
bulge fractions {B/T) of the 15 E+As for which the 
bulge-disk decomposition is rehable (red shaded his- 
togram). The decompositions are not reliable in other 
cases, for example, in those galaxies (EA07, 18, 19) where 
a large portion of the bulge is masked by dust lanes and 
in strongly disturbed systems (EAOlA, EA02, EA12) 
where the distinction between bulge and disk is mean- 
ingless. Except for two disky E+As (EA09 and EA17), 
most (13/15) of the remaining E +As have B/ T's greater 
than the median for Sa's (0.45; iKentI Il985[ ). The me- 
dian B/T of E +As (0.59) is consistent with that of SO 
galaxies (0.63: iKenti 119851 ). E+A galaxies are mostly 
bulge-dominated systems, confirming previous studies 
using bulge-disk decomposition (e.g ., lYang et al] |2004 
IBlake et al .112003 : IB alo eh et al.ll2005h . 



3.3.2. Sersic Profile 

The Sersic i ndex n is often used t o quantify galaxy mor- 
phology (e.g.. iBlanton et al]|2003D because r^/" profiles 
are flexible enough to represent a wide range of profiles, 
from those of disk galaxies (n = 1) to those of spheroidals 
(n = 4) . E+A galaxies in the SDSS have large Sersic in- 
dices (median n = 3), and therefore their l uminosity is 
centrally concentrated (jQuintero et al.|[2004f ). We fit sin- 
gle Sersic profiles to all of our E+A galaxies, and present 
the results in Table 21 This model generally results in 
larger xt values than those obtained with bulge+disk 



decomposition, but both sets of models are statistically 
acceptable. 

In Figure [3f), we plot the distribution of the Sersic in- 
dex n in the R band for 20 E+As, the entire sample 
except for EAOlA, which is too disturbed to fit. Most 
E+As (17/20) have very high Sersic index, n > 5, even 
higher than the n = 4 de Vaucouleurs profile of nor- 
mal ellipticals, indicating that the galaxy's luminosity is 
highly concentrated, presumably due to recent central- 
ized star formation. The exceptions are the three dusty 
galaxies (EA07, EA18, EA19) with n = 1 - 3, where we 
hypothesize that the dust lanes partially mask the bulge 
component. High central concentrations do not necessar- 
ily mean that the galaxies are bulge-dominated systems. 
For example, two of our most disky systems, EA09 and 
EA17, with B/T = 0.17 and 0.35, respectively, have un- 
usually high Sersic indices (n > 10). 

What causes the Sersic index of E+As to be even 
larger than those of early-types (e.g ., [Kelson et al.| [ 



iGraham et 311120011: iGrahamI 120021 )? We consider three 
possibilities: 1) we are being misled by the comparison 
across different studies because the derived Sersic indices 
are strongly dependent on the fitted radial range, the 
choice of fitting method (ID vs. 2D), and the depth of 
the images, 2) the index is dominated by the presence of 
an additional, highly concentrated, central source, and 
3) misidentification of a disk or tidal feature as the shal- 
low wing of the Sersic profile is creating the apparently 
large central concentration. For the remainder of this 
discussion, we presume that the first possibility is not 
the origin of these results, although a more thorough 
cross-comparison among samples is warranted. 

To test the degree to which the central regions of E+As 
affect the Sersic fits, we mask the centers with a 3 pixel 
radius mask and repeat the fitting. In Figure 21;, we 
plot the Sersic indices before and after masking. Most 
(65%) of the galaxies have a smaller n after masking, 
which suggests that a central concentration does play a 
role in causing the high n's. This conclusion is supported 
by positive central residuals in the model subtracted im- 
ages (Fig. [H and H]). EA09 is a special case because 
the bulge is so small that this masking leaves only the 
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disk component. However, Sersic indices in six E+As 
actually increase after the masking. The E+As with a 
larger resulting n after masking generally have ring-like 
substructures in the residual images and negative central 
residuals. Wc conclude that substructures near the cen- 
ter, e.g., bright nuclei, bars, and rings, play an important 
role in generating the highly concentrated light of most 
E+As. 

The last of our three possibilities is that disk or tidal 
features are fitted as the shallow wings of the Sersic pro- 
files. To examine the Sersic index of the bulge component 
alone, we fit a Sersic-l-disk model. In Figure [4l>, we show 
the resulting distribution of Sersic indices for the bulge 
components (red shaded histogram) compared with the 
original, single-component values. The Sersic index of 
the bulge itself is more consistent with that of normal 
bulges (n ~ 3 — 5); however, some bulges still have large 
Sersic indices (n > 5). We conclude that both the cen- 
tral substructures and disk-like features are responsible 
for the unusually high Sersic indices in E-|-As. In any 
case, our E-|-As are more consistent with n ~ 4 than 
n~ 1. 

3.4. Concentration and Asymmetry, Residuals 

So far, we have derived structural parameters by fitting 
smooth symmetric models to the images, even though 
asymmetric features are quite common in the resid- 
ual images (Fig. [T] and H]). To quantify the asym- 
metric features and measure morphologies in a model- 
independent way, we calculate nonparametric measures 
of the galaxy m orphologies: the concentration index C 
(lAbraham et al . 1994) and the rotational as ymmetry in- 
dex A (jSchade et al.iil995l ) . In previous work (jYang et al.l 
I2004D . we u sed the azimuthal F ourier decomposition 
method fe.g.. lRix fc Zaritskvlll995D to investigate asym- 
metric features of the disk-like components, such as the 
lopsideness in two face-on disky E-|-As (EA03 and EA04). 
We use CA indices in this paper, because a large com- 
parison sample for the azimuthal Fourier decomposition 
is not available. 

To determine the C and A indices (collectively re- 
fer red to as CA), we follo w the methodology d e scribe d 
by IConselice etall (|2000[ ) and iBershadv et al] (|2000D . 
First, the Petrosian radius rp is determined from the 
growth curve for circular apertures requiring that ry = 
I{rp)/{I{rp)) = 0.2, where I{r) is the surface brightness 
at the ra dius r and (/(r)) is the mean surface brightness 
within r (jPetrosianI 1 1 976( 1 . We calculate the CA indices 
for image pixels within the "total" aperture defined to 
be twice the Petrosian radius. We measure the asym- 
metry index A by subtracting the galaxy image rotated 
by 180° from the original image and correcting for the 
contribution to the asymmetry arising from background 
noise: 

A.m^n^ El^ol J E l^ol ' 

where / and B represent the image pixels within a cir- 
cular aperture of radius 2rp and the background region, 
respectively. We choose a rotational center that mini- 
mizes A by searching over sub-pixel grids with 0.1 pixel 
resolution. The concentration index C is defined by 



u 



a 
o 
u 




1. 



+ 



0.0 0.1 0.2 0.3 

Asymmetry A 



0.4 



C^51og(^ 
r20 



(3) 



Fig. 5. — CA classification diagram for our 21 E+As as well 
as 113 loc al galaxies with various Hubble types drawn from the 
IFrei et aLl ([1996) catalog. The ellipses, crosses and spirals repre- 
sent early- types (T ^ 0), intermediate spirals, and late- type spirals 
(T ^5), respectively. E+As are represented with filled circles. We 
overlay a dashed line providing a rough division of the C — A plane 
into regions of early- and late-types. E+As populate the CA plane 
widely, as suggested by their morphological diversity. Two-thirds 
of E+As occupy a unique part of the plane: high concentration (C 
> 4.3) and moderately large asymmetry (A > 0.04). These E+As 
have concentrations consistent with those of early type galaxies 
((C) = 4.4 + 0.3) or even higher, as expected from their B/T frac- 
tions and high Sersic indices n. However, they have considerably 
larger asymmetry than ellipticals {{A) = 0.02 ± 0.02) due to dis- 
turbances in their centers, presumably arising from the starburst 
and/or recent galaxy-galaxy interaction. 

where r2o and rso are the radii containing 20% and 80% 
of the total fiux, respectively. Galaxies with a de Vau- 
couleurs profile (n = 4) or a pure exponential disk profile 
have C — 4.50 or 2.79, respectively. 

In Figure [Sj we show the CA classification diagram, 
which is frequently used for t he morphological class ifica- 
tion of high- 2; galaxies (e.g. lAbraham e t al.l [l996[ ). for 
our 21 E+As and for 113 lo cal galaxies with various 
Hubble types drawn from the iFrei et al] (|1996f) catalog. 
The CA indices of the Frei sample were measured in ex- 
actly the same way as those of the E+As. Although the 
CA diagram is not an ideal classification tool (e.g., see 
IConselice et al.l 2000) , it is quite useful in distinguishing 
early from late type galaxies. In Figure [5l intermediate 
spirals occupy a wide space in the CA plane, but early- 
types are segregated from late-types. The CA indices 
have low sensitivity to S /N and spatial resolution, within 
certain limits. For example, Lotz ct al. (2004) show that 
CA indices can be robustly measured within AC < 0.1 
and AA < 0.05 for images with an average (S/N) > 5 
and a s patial resolution of ^ 500 pc per pixel or better 
(see also lConselice et "al]l2000l ). Because our HST images 
satisfy these criteria {{S/N) > 7 within rp and a spatial 
resolution < 100 pc per pixel), a direct comparison be- 
tween the E+As at z '--^ 0.1 with the local Frei sample is 
valid. 

E+As populate a wide range in the CA plane, as al- 
ready suggested by their morphological diversity. We 
overlay a dashed line providing a rough division of the 
CA plane into regions of early- and late-types. Most of 
the dusty E+As (EA07, 12, 13, 18, 19) and the train- 
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wreck E+A (EAOIA) lie below the dividing line and 
have large asymmetry and/or low concentration, presum- 
ably due to extinction. The remaining E+As are located 
above the line, but populate a unique part of the CA 
plane where there is both high concentration C > 4.3 
and moderately large asymmetry A > 0.04. These E+As 
have concentration indices greater than or equal to those 
of early types ((C) = 4.4 ± 0.3), as expected from their 
B /T fractions and high Sersic index n's, but they have a 
considerably larger A than ellipticals {{A) = 0.02 ±0.02). 

Where in radius does the largest contribution to the A 
parameter come from? Faint tidal features do not con- 
tribute significantly to A because 1) the total aperture 
size is not large enough to enclose features at such large 
radii and 2) the asymmetry is normalized using the to- 
tal galaxy luminosity (eq. so little weight is placed 
on low surface brightness tidal features. In the local 
galaxy sample, small scale features such as spiral arms 
an d flocculence in the dis k are the main contributors to 
A (jConselice et al.l l2000fl . However, such features are 
not visible in the disks of our E-t-As (except for in EA03 
and EA19), which are sm ooth (see also the NIR imaging 
from lBalogh et aLll2005f l. Therefore, the large asymme- 
try of E-I-As is probably due to relic structures within a 
few effective radii arising from the starb urst and/or re - 
cent galaxy-galaxy interaction. Yamauc hi et al.l (120051 ) 
reached a similar conclusion from their analysis of 22 
SDSS E-hAs. 

To further investigate the nature of the asymmetry, 
we compare the asymmetry profiles, i.e., the asymme- 
try index interior to radius r, A{< r), of E-|-As and 
late-type spirals (T ^ 5). Except for three E-l-A galax- 
ies (EAOIB, 03, 11), E-I-As have decreasing or flat A{r) 
profiles. Conversely, the late-type galaxies tend to have 
increasing asymmetry profiles, because spiral arms and 
flocculence in the disk are increasingly included as the 
aperture increases. We conclude that the asymmetry in 
E-|-As originates mostly in the central parts (within a few 
Te), presumably due to unrelaxed structures or residual 
dust from the starburst and/or recent merger. These 
E-t-As would be classified as "normal" elliptical galaxies 
or early-type spirals if they were observed at high red- 
shift, unless the disturbances in their innermost regions 
(within a few kpc) or the low surface brightness features 
with /J, ~ 24 — 25 could be clearly identified. In Table [31 
we list the aperture sizes (r^), CA indices, and relevant 
aperture magnitudes. 

4. COLOR PROFILES 

The radial color profile of a galaxy depends on its dust 
content and the spatial distributions, ages, and metallic- 
ities of its stellar populations, which in turn depend on 
the evolutionary history of the galaxy. For example, if 
galaxy-galaxy interactions are responsible for an E+A's 
recent starburst, the young stellar population is expected 
to be centrally concen trated (e.g., Mihos & Hcrnquist 
119941 : iBekki et al.ll2005l ) and to produce a positive color 
gradient, i.e., a redder color with increasing radius. On 
the other hand, if the E-l-A arises from a truncation of 
continuous star formation, due to a mechanism such as 
ram pressure stripping in the cluster environment, then 
it may have a uniform color profile (see the exten sive 
discussions in lCaldwell et all 1 19991 : [Rose et al.ll2001h . 

Current observational evidence supports the model 



of centralized star formatio n in E-l-A ga l axies. Us- 
ing long slit spectroscopy, [Norton et al.l (|2001f) find 
that young stars are more centrally concentrated than 
older ones in our E-|-As. Even for E-|-As in clusters, 
[Bartholomew et aTl (2001) show that on average they be- 
come slightly bluer towar d the center than norma l early- 
type galaxies. Recently, [Yamauchi et al.[ ([2005f ) found 
that a substantial fraction of SDSS E-|-As have such a 
positive color gradient. However, these studies are some- 
what limited due to the galaxies' small angular extent 
and/or the effects of ground-based seeing. In this section, 
we present high resolution internal color distributions of 
our E-|-As and show that indeed a significant fraction 
have positive color gradients and sometimes distinct blue 
cores. Then, we investigate whether the color profiles of 
E-|-As can evolve into those of early- type galaxies. 

4.1. Positive Color Cradients in E+A Galaxies 

We show the (-B435— i?g25) radial color profiles of 20 
E-|-As in Figure [H excluding EAOlA because it is so 
irregular. To enable direct comparison between the 
profiles, the (-B439— i?702) colors were transformed into 
(^435— -^625), and all colors were X-corrected with re- 
spect to 2; = 0.1 to minimize the X-correction uncertain- 
ties. We refer to these colors as (B — i?)°-^. Note that 
we do not apply different K-corrections along the color 
gradients because these corrections are small (typically 
0.02 mag) compared to the overall K-correction values. 
We measure the radial color profiles using the surface 
brightness profiles obtained with the ELLIPSE task. Be- 
fore fitting ellipses to the isophotes, we correct for the 
different point spread functions in the B and R bands 
by convolving the sharper R band images to match the 
B band images. The smoothing kernel is determined us- 
ing unsaturated stars in both images. We extract surface 
brightness profiles in the R band as described in i )3.3[ and 
then extract the B band profiles using the same isophotes 
as in the R band analysis. 

Note that this PSF-matching is critical in order to in- 
vestigate the innermost color profiles of E-l-A galaxies 
and to identify color structures such as blue cores that 
can be as small as a kpc (see also [Menanteau et al.l ([2004[ ) 
for a comprehensive treatment of PSF issues). Figure [3 
shows the effect of this PSF correction. Assuming an 
intrinsically uniform color distribution, i.e., a zero color 
gradient, we generate model B and R images with the 
best fitting GALFIT models for our 20 E-fAs, and then 
convolve these with the appropriate B and R PSF's to 
generate the "observed" B and R images. In the upper 
panel, we plot the observed color profiles when we do not 
match PSFs. We plot the measured color profiles after 
PSF-matching in the lower panel. The improvement is 
obvious and the measured color profiles deviate from the 
intrinsic flat proflle by no more than ~ 0.02 magnitudes. 
The extracted profiles tend to be biased toward redder 
cores, but the effect is negligible. 

We find that all E+As have globally blue colors, ex- 
cept within dusty regions, and so conclude that recent 
star formation extends acros s the face of these galaxies 
(consistent with the results of [Franx[[1993[ : [Caldwell et al.l 
1996; Nort on et al.l 2001). In Figure [6l we overlay the 
model {B — i?)°-^ colors of 5 to 10 Gyr-old simple stellar 
populations which are consistent with those of early type 
galaxies for comparison. 
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Fig. 6. — Redshifted radial (B435— i?625)'''^ color profiles of the 20 E+As. All the colors are K-corrected to z = 0.1. The color profiles are 
as diverse as the morphologies of E+As, ranging from negative color gradients (EA02, 04, 07, 12, 19) to positive color gradients including 
compact blue cores (EA05, 06, 09, 14, 16, 17). The large fraction (70%) of positive color gradients suggests that many E+As had a central 
burst of star formation and therefore that the E+A phase is triggered by galaxy interactions/mergers that funnel gas toward the center 
of the remnant. The shaded region indicates the (B435— i?g25)'^'^ colors of a 5—10 Gyr old SSP with metallicity Z = 0.004 {1/5Zq). The 
color of an SSP with solar metallicity is redder and varies from 1.86 to 2.04 between 5 and 10 Gyr. Recent star formation is not limited to 
the central region, as evidenced by the fact that the color is bluer than the shaded region at most radii. 
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Fig. 7. — Effect of the PSF correction. (Top) Measured color 
profiles for model galaxies witli a uniform color distribution, i.e., a 
zero color gradient, without convolving the PSF-matching kernels 
with the R band images. The non-zero color gradients arising due 
to differences in the PSF function are clear over the entire range 
of radii depending on the model galaxy. (Bottom) Color profiles 
extracted after the PSF-matching. The measured color profiles 
deviate from the fiat profiles by no more than ~ 0.02 magnitudes. 
Our color gradient measurement is valid even in the innermost 
region. 

While E-|-As are relatively blue overall, their inter- 
nal color variations are surprisingly diverse, as suggested 
by their morphological diversity (see Table (HI). Twelve 
E-|-As (57%) have positive color gradients, i.e., they be- 
come redder with increasing radius, while five E-|-As 
(20%; EA02, 04, 07, 12, 19) have negative color gradi- 
ents. The remaining five E-|-As have relatively flat color 
profiles or a mixture of positive and negative color gra- 
dients. Half of the E+As with positive gradients exhibit 
blue cores, which we discuss below. This diversity echoes 
the variation in US abso rption line strengt h profiles ob- 
served in cluster E+As (jPracv et al.ll2005f) . We list the 
overall color morphologies of E-l-A galaxies in Table [6] 
within broad categories: positive (including blue core), 
negative, or flat color gradient. 

Among the E-f As with negative color gradients, three 
(EA07, 12, 19) show clear dust signatures, such as dust 
lanes and irregular filamentary structures in the two- 
color composite images. The clear detection of dust fea- 
tures in the ACS' sample is mainly due to the excellent B 
band sensitivity, and so we cannot exclude the possibility 
that similar dust exists in the WFPC2 sample galaxies 
with negative color gradients (EA02 and 04). Therefore, 
although there are other possible explanations, including 
a trun cated IMF during the starburst (see iBekki et all 
(|2005f )). we speculate that the red cores in these galaxies 
arise from increasing dust extinction toward the center. 

To quantify the color gradients, A{B — i?)/Alogr, we 
perform linear least-square fits to the radial color pro- 
files. Because of the highly variable color profiles and 
compact blue cores in some E-|-As, a single color gradi- 
ent is obviously not the best representation. Therefore, 
we measure the color gradients with two piecewise linear 
fits, allowing the break radius i?brcak between the inner 
and outer slope to change freely. Our only constraint 
is i?brcak ^ 0.3 kpc (~ 4 pixcls) to prevent measuring 
break radii that might be seriously affected by the PSF. 
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Fig. 8. — Distribution of the E-l-A color gradients A(B — 
-R)/Alogr. The solid and shaded histograms represent the gra- 
dients measured from the two-line fits and single line fits, respec- 
tively. Most E+A galaxies have positive color gradients, i.e., the 
color becomes redder as the radius increases. Note that there are 
some very steep [A(iJ — iJ)/Alogr > 0.6] color gradients due to 
the blue cores. The dashed curve represen ts the distribution of 
the color gradients of local E/SOs llWise &: S ilva 1996) with arbi- 
trary normalization. E/SOs have a narrow range of color gradients 
[A(B - R)/Alogr = -0.09 ± 0.06]. These color gradients are in- 
terpreted as being the result of metallicity gradients. 

In Table [6l we list the color gradients and break radii 
from the piecewise linear fits, as well as the color gradi- 
ent derived from the single line fit. The color gradients 
are measured in (i?435— -R625) at z = 0.1 to minimize the 
systematic errors that may arise from the K-corrections 
and color transformations. For EAOl-05, the WFPC2 
sample, we measure the color gradient in (i?439— i?702) 
at z = 0.1 and convert it to (i?435— i?625) at z = 0.1 
using the relation A(B435-i?625) — ci A(i3439-i?702) in 
Table [21 We find that the broken linear fits significantly 
improve the fits (x^s) in aU but three E+As (EAOIB, 03, 

The distribution of color gradients is shown in Figure[8l 
where the solid and hatched histograms indicate the color 
gradients for single and double line fits, respectively. The 
range of color gradie nts fo r early-type g alaxies (E and 
SOs from lFranx et al.lll989HPeletier et al.lll990( ) is shown 
for comparison (dashed curve). The majority of E-t-As 
(~ 70%) have a positive color gradient that deviates from 
the slightly negative profile typical of early-type galaxies. 
We conclude that the young stellar populations are more 
concentrated than the underlying old populations. This 
conclusion is independent of whether the single or double 
line fits are used and is qualitatively consistent with that 
of .Yamauchi et al.. (2005). 

4.2. Blue Cores 

The most interesting morphological feature in certain 
E+As is the compact, almost stellar-like, blue core that 
is distinct from the other parts of the galaxy. We de- 
fine blue core galaxies as those that have a color vari- 
ation, 5{B — R), within the inner 2 kpc region that is 
greater than 0.3 magnitudes. This defini tion is compara- 
ble to the original criteria adopted by JVI enanteau et al.l 
(|200H ). S(y - /)obs > 0.2, for morphologically selected 
early-type galaxies. Six E+As (EA05, 06, 09, 14, 16, 17) 
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Fig. 9. — Examples of color gradient evolution for two E+As with positive gradients (EA3 and 10) and for two blue-core E+As (EA16 
and 17). The solid and dashed lines represent the current color profiles and the subsequent evolution, respectively. The ranges of color 
profiles after 5 and 10 Gyr are represented with solid and hatched regions, respectively. We adopt Zyoung = 0.02 for EA03 and EA16 (left), 
■ 0.05 for EAIO and EA17 (right). The chosen burst ages (tburst) can produce the whole range of colors in each galaxy. The 



and Zy 

positive color gradients in E+As can be inverted within 
with that observed for E/SOs. 

satisfy our definition and have characteristic core sizes 
of i?bioak = 0.4 - 1.4 kpc. Within these break radii, 
the color gradients of blue core galaxies are very steep, 
A{B — R)/ Alog r> 0.6. The slopes become relatively flat 
at radii larger than the break radii. Note that we clas- 
sified EAOIB as one of the blue core galaxi es based on 
its \V FPC2 {B — R) colors in previous work (|Yang et all 
|2006| ): however, EAOIB is not counted as a blue core 
galaxy in this paper because its color profile does not 
have a break radius. 

Though the origin of blue cores is not fully understood, 
they are common in early-type galaxies at higher red- 
shifts (z > 0.5), when field spheroids are presumably still 
assembling. For example, 30% of the morphologically- 
selected elliptical galaxies in the Hubble Deep Field 
North h ave color inhomogene iti es, mostly d ue to blue 
cores [M enanteau et al.ll2001h . iTreu et all (|2005D find 
that ^ 8% (14/165) of early-type spheroidals in the Great 
Observatories Origins Deep Sur vey North (GOODS-N) 
have blue cores. I Lee et aP ()200 6f) also find "blue clumps" 
in almost half of the early-type galaxies with blue colors 
in the GOODS North and South fields. The blue-core 
E-|-As in our sample have relaxed morphologies and bulge 
fractions of 0.2 — 0.6. Given that the CA classification 
scheme identifies these galaxies as early type galaxies, 
the blue-core E-|-As would be classified as early-types if 
they were observed at high redshift. Therefore, these 
blue core E-l-A galaxies might be the local analog of the 
blue core spheroids often found at high redshift. The 
importance of blue core E-|-As for better understand- 
ing galaxy-galaxy mergers and the origin of the black 
hole mass-bu lge velocity dispersion relation for galax 
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5 Gyr, evolving into a range [—0.05 < A(_B — _R)/Alog r < —0.16] consistent 



I2000f) is described bv lYang et all (|2006f) . 

4.3. Evolution of Color Gradient 

It is now well established that E/SOs in the local Uni- 
verse have negative color gradients [A{B — i?)/Alogr = 
—0.09 ±0.06 mag dex~^ in Fig. [5], which originate from 
their metallicity gradients; their stellar populations be- 
come more metal ric h and redder toward the center (e.g., 
iTamura et ani2000[ ). If E-hAs evolve into E/SOs, how 
can their diverse color morphologies converge into such 
a narrow range of color gradients within a few Gyr? In 
particular, is it possible that the observed positive color 
gradients invert into the negative color gradients typical 
of early-types? 

For those E-|-As whose negative color gradients arise 
from dust, it is difhcult to determine the underlying color 
gradient s without incorpora ting radiative transfer calcu- 
lations (|Wise fc Silvalll996[) . Nor do we know how the 
dust content will evolve. However, were the dust in these 
E+As to disperse somehow, the negative color gradients 
[A{B — R)/A\ogr < —0.3] would flatten and be more 
consistent with those of E/SOs. 

For those E+As with positive color gradients — which 
are most of the sample — the centralized starburst will 
naturally lead to metallicity and color gradients if the 
young stellar populations are more metal rich than the 
underlying old populations. Although violent mixing 
during the merger could dilute any previously established 
metallicity gradients, numerical simulations suggest that 
metallicity gradient s might be regenerated by central- 
ized s tar formation (jMihos fc Hernqui"stlll994l : iKobavashU 
[200l . 

To investigate the evolution of the color gradients of 
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E+As quantitatively, we take the current color gradi- 
ents and evolve the stellar populations using simple as- 
sumptions. We assume that the current color gradients 
are due solely to burst-strength gradients and that the 
young populations are coeval, i.e., that the last star- 
burst was instantaneous. We adopt a uniform metal- 
hcity (Zoid = 0.004) for the old population, and consider 
^young — 0.02 and 0.05 for the young population. Then, 
for each galaxy, we choose a single post-burst age (the 
time since the starburst ended) that produces the whole 
range of colors observed in that galaxy, and derive the 
burst strength profile that reproduces the color profile. 
As older post-burst ages are selected, the burst strength 
gradient increases for a given color, and therefore larger 
metallicity gradients are introduced. We evolve these 
stellar populations passively for 5 and 10 Gyr ignoring 
further dynamical evolution, and measure the resulting 
color gradients. 

In Figure [HI we show the color gradient evolution for 
two E-t-As with positive gradients (EA03 and 10) and 
for two with blue-cores (EA16 and 17). The solid and 
the dashed lines represent the current color profiles and 
their subsequent evolution, respectively. Taking into ac- 
count measurement errors, we show the range of color 
profiles after 5 and 10 Gyr with the solid and hatched 
regions, respectively. We adopt .^young — 0.02 for EA03 
and EA16 (left panels), and Zyoung — 0.05 for EAlO and 
EA17 (right panels). 

In these examples, the positive color gradients can in- 
vert within ~ 5 Gyr and evolve to a color gradient range 
[-0.05 < A(B-i?)/Alogr < -0.16] consistent with that 
of E/SOs. Larger metallicity differences between the old 
and young stellar populations and older post-burst ages 
will result in steeper metallicity and color profiles after 
5-10 Gyr. While it is hard to predict the true color gra- 
dient evolution of the E-|-As without exact knowledge 
of the post-burst ages and metallicities, our calculation 
shows that the current color profiles, most likely arising 
from a centralized starburst during a galaxy-galaxy tidal 
interaction or merger, can evolve into those of early-type 
galaxies. 

5. E/SO SCALING RELATIONSHIPS AND E+AS 

In previous sections, we have shown that some as- 
pects of E-l-A morphology [B/T, concentration C, Sersic 
index n) are consistent with those of current early- 
type galaxies. We have also presented models or ar- 
guments in which other aspects now inconsistent with 
current early-type properties (blue cores, tidal features, 
and positive color gradients) could, via evolution, be- 
come consistent with the properties of current early- 
type galaxies. To further investigate whether E-l-A 
galaxies truly evolve into present-day early-types after 
a few Gyr, when the discrepant morphological features 
might disappear, we test whether E-|-As will eventu- 
ally lie on the various scalin g relationships of e arly- 
type galaxies. For examp le, iTacconi et al.l (|200^ and 
iRothberg fc Joseph! (|2006l ) show that late-stage ultra- 
luminous infrared galaxy (ULIRG) mergers and optically 
selected merger remnants (other possible progenitors of 
present-day E/SOs) are located very close to or on the 
K-band Fundamental Plane (FP) of early type galaxies. 
We start with the simplest forms of the scaling relations, 
i.e., the Faber- Jackson relation (FJR; iFaber fc Jackson! 



I1976D and the Kormendy relation (KR: lKormendvlll977f) . 

5.1. Faher- Jackson and Kormendy Relations 

In Figure (TOl we show the Faber- Jackson relation 
[L — a'^) and the Kormendy relation (re — (/ie)) of E/SOs 
in comparison to those of E-|- As. We draw the E/SO 
comparison sample from lJorgensen et al.l (|1996f ). because 
their photometric bands (Gunn r and B band) closely 
match our HST WFPC2 and WFC filters. For the 
E+As, we adopt th e velocity dispersion s of the old com- 
ponents (uoid) from lNorton et all (|200lD . To disentangle 
the kinematics of th e young and old stellar populations, 
iNorton et al.l (|2001h simultaneously fit the velocity pro- 
file and the relative contributions of different stellar com- 
ponents using longslit (1.5" x 6. 6"aperture) spectra. For 
three E-|-As (EA06, EA17, and EA18), we have only up- 
per limits on CToU, because a large fraction of the light is 
contributed by a young A-type population {Ja in NOl) 
and thus the velocity dispersion of the old K-type pop- 
ulation is not measured reliably. Although we include 
these galaxies in Figures [10! and [121 we exclude them 
from the following analysis. For the structural param- 
eters, we adopt r^., (/ie) and the total magnitude Mtot 
from our de Vaucouleurs fit (Figure [TT]). While a single 
de Vaucouleurs profile is not an ideal fit/model for com- 
plex morphologies (e.g., EA02) and/or for disky systems 
(EA09 and EA17), we justify the choice of rg from the 
r^/"* fits as a simple measurement of the half light radius 
of the galaxy. Indeed, we find that the luminosity within 
rp given by the GALFIT fit agrees with a circular 
aperture measurement of the half light radius to within 
10-20%. 

We reproduce iNorton et al.l (|2001f )'s result that there 
is a strong correlation between E-l-A magnitude and ve- 
locity dispersion, and find that it is roughly parallel to 
the E/SO Faber- Jackson relation with some offset (left 
panels in Figure fTOj) . This relation implies that E-l-A 
galaxies already harbor dynamically relaxed old popu- 
lations an d that they are pressure-supported systems. 
While Nor ton et al.l |2001h had to perform several steps 
to transform the LCRS Gunn/Kron-Cousins magnitude 
into Gunn r magnitude to compare the properties of 
E+As with those of E/SOs, we determine more directly 
the relative offsets between the E-|-As and E/SOs with 
our HST photometry. We measure the offset by fitting 
the L oc (T^ (i.e., Mtot oc 10 log u) relation to the E+As 
and the J96 sample. 

On average, E-l-A galaxies are 1.31 (±0.73) and 1.64 
(±0.71) mag brighter than the E/SOs in the Gunn r 
and B bands, respectively. This offset in Gunn r is 
almost twice as lar ge (~ 0.6 mag) as that derived by 
INorton et all (|2001h . This difference arises not only be- 
cause we use improved photometry and magnitude trans- 
formations, but also because the functional forms used 
to calculate the offset are different. 

These offsets can be interpreted as either that E+As 
have an elevated luminosity for a given velocity disper- 
sion or that they have an unusually low velocity disper- 
sion at a given luminosity. To investigate which interpre- 
tation is more correct, we compare E+As to E /SOs on the 
purely photometric scaling relation, i.e., the iKormendvl 
(|1977f ) relation (hereafter KR), between the effective ra- 
dius (re) and the mean surface brightness within r^. The 
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Fig. 10. — Faber-Jackson and Kormendy relations for E-f A galaxies in the B {bottom) and Gunn r bands (top). The small squares 
and large circles represent cluster E/SOs in the J96 sample and our E+As, respectively. The large filled squares are the E+As for which 
we have only an upper limit on the velocity dispersion of their old stellar population, a^id from NOl. In the left panels, the solid and 
dot-dashed lines represent the best-fit Faber-Jackson relations (L oc ct*) for E-I-As and E/SOs, respectively. On average, E-l-A galaxies are 
1.31 (±0.73) and 1.64 (±0.71) mag brighter than E/SOs in the Gunn r and B bands. In the right panels, the solid and dot-dashed lines are 
the Kormendy relations with fixed slope for E+As and E/SOs, respectively. The average offsets between E+As and E/SOs are 0.98 ± 0.45 
and 1.57 ± 0.52 mag in r and B, respectively. The dashed line is the boundary set by the limiting magnitude. The dotted line is the 
sharp p hysical boundary tha t divides the plane into the E/SOs and the so-called exclusion zone where dynamically relaxed systems are not 
allowed ^Bender et al.|[l99l) . 

conclude that elevated luminosities, probably due to the 
recent star formation, are responsible for the offsets of 
the E+A scaling relations, rather than unusually low ve- 
locity dispersions. While the FJR and KR suggest that 
roughly one magnitude of fading may bring E+As down 
to the E/SO scaling relations, the large scatter makes 
such an analysis difficult. 

5.2. Fundamental Plane 

The fundamental plane, hereafter FP, is an empirical 
scaling relation between the effective (or half-light) ra- 
dius Te, the central velocity dispersion cr, and the mean 
surface brightness (/e) within re for earl y type galaxies 
(jPiorgovski fc Davislfl987l : lDressleilll987| ). Although the 
physical origin of the FP is not fully understood, both 
cluster and field early-type galaxies follow this relation 
with remarkably small scatter (e.g., ~ 0.1 dex in log re in 
the Gunn r band; J96), and an extension of this formal- 
ism appears to fit all spheroids ranging from Galactic 



right panels in Figure [TUl show the KRs for the E+As and 
E/ SOs, the latter of which is corrected from the original 
KRiKormendv (19 7^) by assu ming an average color, B—r 
~ 1.20 (|Fukugita et al.ll 19951) . The slope of the relation 
is known to be insens i tive t o the band. For example, 
iHamabe fc Kormendvl (|1987f ) find a slope of 2.98 in the 
V band compared to 3.02 in the B band. 

In the KR plane, most of E+As, except for a few disky 
ones, have a higher surface brightness than E/SOs at a 
given effective radius. The average offset between E+As 
and E/SOs is 0.98 ±0.45 and 1.57 ±0.52 mag in r and 
B, respectively. The offsets appear to increase as the 
galaxies become smaller; roughly half of the E±As do 
not overlap with the B and r band E/SO loci. The lower 
boundary (dashed line) of the E/SO region is set by the 
limiting magnitude of ~ —20.45. The upper bound- 
ary is not due to selection effects (see ^5.3p . Because 
the surface brightness offsets in the KR are consistent 
with the magnitude offsets obtained from the FJR, we 
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The structural 



dSph galaxies t o the intracluster st ellar component of 
galaxy clusters (jZaritskv et al.|[2006l ). 

To examine the FP of E+As and compare it with that 
of E/SOs, we adopt th e definiti on and methodolog ies of 
iJorgensen eFall (|1996D (see also lKelson et al.|[2000l ). The 
fundamental plane is defined as 



log Te = a log cr + /3 log (/e) + 7 , 



(4) 



where Vg, a, and (/e) are the effective or half-light radius 
in kpc, central velocity dispersion in km s~^, and mean 
surface brightness within in units of Lqpc^^. For the 
E-l-A structural parameters, we adopt and (fig) from 
the single de Vaucouleurs fits to our HST images (Fig- 
ure [TT|) . We transform the mean surface brightness in 
our HST bands into Gunn r and Johnson B using our 
transformation relations (eq[T]) and the coefficients in Ta- 
ble [21 We convert to units of L0pc~^ using the relation 
(/g) = — 0.4((/ie) — constant), where the constants are 
26.4 and 27.0 for Gunn r and Johnson B magnitudes, 
respectively. We do not apply a color-dependent correc- 
tion to Tg, because the E/SO FP is only weakly s ensitive 
to color (Pahre et al. 1998; Bernardi et al.ll2003^ . 

5.3. Location of E+As on the FP 

We show the location of our E-f As on various pro- 
jections of the J96 FP in the Gunn r and B bands in 
Figure [121 Figures [T2b and [T2b present the face-on view 
of the E/SO FP in r and B, respectively. The dashed 
line is the boundary set by the limiting magnitude of 
Mr ~ —20.45, which corresponds to the constant mag- 
nitude line (dashed line) infTOb andfTOli. The boundary 
set by the dotted line divides the plane into the E/SO FP 
and the so-called zone of exclusion ( ZOE), where dynam - 
ically relaxed systems do not reside ([Bender et al.lfl99^ . 
The physical origin of this boundary is not fully under- 
stood, although it might correspond to the upper limit 



of the velocity disper sion function of early-type galax- 
ies (cr < 350 km s-i; iTreu et all [20061) . In the face-on 
projection, E-|-As span the same region defined by the 
normal early-type galaxies in the r band; however, half 
of our E-|-As violate the zone of exclusion in the B band. 
We plot vectors parallel to the FP parameters to find out 
which parameters drive this shift in the r band. E+As 
are shifted from the E/SOs along the direction of increas- 
ing (/g), as we concluded earlier from the Kormendy rela- 
tion. Given that E-|-As already populate the same region 
as E/SOs in the r band and that the (/e) offset is larger 
in B, this violation of the ZOE must be due to the tem- 
porary increase of the surface brightness caused by the 
recent starburst. 

Figures [T2b and[T2lj present the edge-on projection of 
the E/SO FP, defined as x — log re and y = alogcr -I- 
/31og(/e), in r and B, respectively. In this projection, 
E-|-As are distinct from E/SOs. The existence of a tight 
correlation among E+As indicates that they follow their 
own scaling relation. The large offset between the E/SOs 
and E+As suggests that the E+As' current stellar con- 
tent is different. Although the mechanisms that trigger 
the E+A phase and stop star formation are not fully un- 
derstood, any model must be able to explain the relative 
tilt and offset of the E+A FP. In the following section, 
we discuss how the observable properties are related to 
physical ones, such as mass-to-light ratio M/L, and pos- 
sible interpretations of the tilt and offset. 

5.3.1. FP Offsets: Evolution of the FP Zero Points 

E+A galaxies stand apart from the E/SO FP in the 
edge-on projection (Figures [T2b and[T2ti). The difficulty 
in interpreting this difference is that any, or all, of the 
three plotted parameters or the implicit parameter, M/L, 
may be evolving as the E+A evolves. The simplest op- 
tion is to assume that only M/L evolves. In this case, 
the offset between the two populations can be directly 
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Fig. 12. — Face-on and edge-on projections of the fundamental plane in the B {bottom) and r {top) bands. The small squares and 
large circles represent the E/SOs in the J96 sample and our E-I-As, respectively. The large filled squares are the three E+As (EA06, EA17, 
and EA18) that have only upper limits on Cold from NOl. Of all the sample galaxies, these three E-I-As have the largest light fraction 
contributed by the young A-type population {Ja in NOl). {a,){c) Face-on projection of the FP. The dashed line is the boundary set by the 
limiting magnitude. The dotted line is the sharp physical boundary that div ides the plane into the E/SO FP and the so-called exclusion 
zone where dynamically relaxed systems are not allowed IIBender et al.|[l992l) . {b){d) Edge-on view of the FP. The solid and dot-dashed 
lines represent the best-fit planes for E+As with and without fixed slopes (a, respectively. E-I-As follow their own scaling relations and 
have significant offsets from the E/SOs. 



related to a difference in M/L 

M _ 



Alog- 



A7 



(5) 



where the FP offset A7 is measured from the difference 
between the intercepts of two parallel lines in Figures 
[H) andlH^. We find A7 = -0.47 ± 0.12 in the Gunn 
r band, which implies Alog (M/L) = -0.57 ±0.14. In 
other words, E+As have, on average, a M/L that is 3.8 
times smaller than that of E/SOs. Note that this aver- 
age offset has a smaller error than the offsets measured 
from the FJR and KR. However, as we have shown previ- 
ously in 21 the structural properties are affected by the 
E+A phase (the color gradients are altered, some E+As 
have blue cores) and so assuming that only M/L evolves 
cannot be completely correct. Our expectation is that 
Te will be smaller in E+As because of the centrally con- 
centrated star formation, and hence that Ve will increase 
(and (/e) decrease) as the E+A evolves, leading to a com- 
plex evolution on the FP. Therefore, the M/L difference 



discussed above is the minimum amount of fading that 
E+As should experience to settle onto the E/SO locus. 
By examining the tilt of the FP, and other projections 
of the scaling laws, we explore whether we can proceed 
beyond this simple model and the potential complexity 
of evolution on the FP. 

5.3.2. FP Tilts: Variation of M/L along the FP 

Given the range of possible evolutionary paths on the 
FP, it is rather surprising that there is a tight relation 
for E+As on the E/SO FP edge-on projection. Fur- 
thermore, the difference appears to be both a relative 
shift and a tilt. To measure the relative tilt of E+As 
with respect to E/SOs, we find the best-fit plane in 
the form of equation [J] through the E+As using an or- 
thogonal fitting method that minimizes the normal dis- 
tances to the plane. We adopt this fitting method to 
enable direct comparison with the fits from J96. J96 
find a = 1.24+0.07, f3 = -0.82+0.02 for Gunn r and 
a = 1.20+0.06, f3 = -0.83+0.02 for Johnson B using 
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226 and 91 cluster galaxies, respectively. Using 16 E+A 
galaxies^ and taking an equal weight for each galaxy, we 
find that the best-fit FPs are 



^1.13±0.10^ 
^1.09±0.08/, 



-0.62±0.07 



-0.59±0.06 



-0.08ifc0.10 



-0.11±0.07 



(6) 
(7) 



for the Gunn r and Johnson B bands, respectively. The 
errors in a, /3, and 7 are estimated using the bootstrap 
method. The rms scatter around these planes is 0.10 and 
0.10 in logTe for Gunn r and Johnson B, respectively. 
The parameters (a, /3, 7) for Gunn r and Johnson B 
agree within the uncertainties, because the FP param- 
eters are generally not highly sensitive to the bandpass 
(J96). Therefore, we adopt the (a, /3, 7) derived us- 
ing the Gunn r data in the following discussion. While 
a = 1.13 for the E-fA FP agrees to within the errors 
with that of the E/SO FP [a ^ 1.24) (the E-hA FP 
is still viewed edge-on in the projection in Figure [T2b ) . 
(3 = —0.62 for the E-|-As is ~ Scr discrepant from that 
of the E/SOs (/3 = —0.82). In the context of interpreting 
differences between the FPs as primarily due to M/L, the 
tilt implies a mass dependence. A similar tilt of the FP 
of merger remna nts is observed relative t o the K-band 
E/SO FP (Fig. 1; iRothberg fc JoseDhll2006h . 

The tight correlation within the E-|-A FP suggests that 
M/L is correlated with the structural parameters. The 
scatter around the plane appears to be comparable to the 
typical scatter of E/SO FPs (~ 0.1 dex). Does this mean 
that E-|-As, many of which are clear merger remnants, 
settle into a relaxed state on a short timescale? Or is it 
possible that the correlation is the result of a selection 
effect? Our selection of galaxies with strong H(5 absorp- 
tion tends to choose galaxies with stronger bursts and/or 
the younger ages; therefore it is possible that our sample 
represents only a certain part of the FP, i.e., near the 
exclusion zone. To address this issue, a larger sample of 
E-|-As covering a wider range of strengths is required. 

Using the best fit FP parameters, M/L for E+As can 
be expressed as 

M/L cx M•'■V^■^ (8) 
compared with M/L cx Af " 24^0.02 f^j. e/SOs. Unlike in 
the E/SO FP, where M/L is a function mostly of mass, 
M/L for E-|-As depends on the effective radius as well 
as the mass of the system. E-|-As with smaller effec- 
tive radii have smaller M/L, possibly resulting from a 
stronger burst and/or shorter time elapsed since the star- 
burst. Because M/L is proportional to {re{Ie))~^ at a 
fixed velocity dispersion, one might expect that a smaller 
Te would lead to a larger M/L, but the observed trend is 
the opposite. Therefore, the surface brightness of E-|-As 
must increase faster than the mass or size of the galaxy 
decreases. This inference is consistent with the observed 
KR (Figure [TO|l. in which (/ie) deviates further from the 
E/SO relation for E-|-As with small effective radii. 

If M/L does indeed vary among E-|-As, then one might 
expect a correlation between the deviation of M/L from 
E/SOs and the color of the galaxy. In Figure [T3l we plot 
AM/L against {B — R) color, and against other param- 
eters (mass, Te, Sersic n, and CTom)- We find a strong 
correlation between AM/L and color, but no significant 

2 This excludes EAOIAB, 07, 17 and 18 which do not have a 
measured cr^i^. 
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Fig. 13. — Correlations between the deviation of M/L from 
E/SOs for individual E+As and the {B — R) color, dynamical mass, 
Sersic index n, velocity dispersion of the old stellar population 
cToitj, and the effective radius re. The squares and circles represent 
parameters from the B and R bands, respectively. AM/L is tightly 
correlated with {B — R) color, but not with n or a. We conclude 
that the tilt, i.e. the variation of M/L along the E-fcA FP, is mostly 
due to stellar population variations. E-f As with smaller effective 
radii have younger stellar populations, on average, and therefore 
smaller M/L. 



correlation with n 01 a. A Spearman rank test shows 
that only {B — R) color correlates with AM/L at a sig- 
nificance level higher than 99.5% for both the r and B 
bands. We conclude that the tilt, i.e., the variation of 
M/L within the E-l-A FP, does reflect a real variation 
of the stellar population. E-f As with smaller effective 
radii have younger stellar populations, on average, and 
therefore, smaller M/L. 

The remaining question is what does the small re phys- 
ically represent? Does it indicate that the intrinsic size 
of the galaxy is small or that there was a strongly con- 
centrated starburst? We compare the effective radii in 
B and r to address these questions. In blue-core E-f As 
or E-fAs with the positive color gradients (~ 75% of 
E-f As), the effective radius tends to be smaller in the B 
band (r^B) than in the r band (rgp;). In these galaxies, it 
is likely that will eventually become close to r^n as 
the light from the young stars in the central region fades 
and the color gradients ffatten. However, most galaxies 
have fairly small differences, |logreB — logreRj < 0.30, 
which is clearly not enough to produce the tilt observed 
over two orders of magnitude in effective radius. 

To further test whether intrinsically small galaxies 
have smaller M/L's, we repeat the orthogonal fitting, but 
this time exclude galaxies with jlogreB — logreR] > 0.15, 
i.e., for which the two effective radii differ by more than 
25%.^ For this subsample, the color gradients are mod- 
erately flat and there is little dust, so strong evolution 
in the effective radius due to stellar evolution is unlikely. 
The fitting of this subsample results in the same FP slope 
as for the full sample, implying that intrinsically small 
and therefore less massive galaxies have smaller M/L ra- 
tios. 

Blue-core galaxies analogous to our blue core E-f As 
(§4.2) are found among early-type galaxies at high-z 

^ We boost the statistics of this test by combining the B and 
R band sample assuming the same FPs parameters for both bands 



(see i|5.3.2l l. 
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(jvan Dokkum fc Ellis! [200l iTreu et al.ll2005f) . and their 
rest-frame B-band M/L is estimated to be much smaller 
(roughly 4x) than that of normal E/SO galaxies at the 
same epoch. As a cautionary note, we examine how 
the difference between the blue and red affects the 
measurement of M/L , especially for blue core galaxies. 
When we use the blue to derive M/ L for our blue core 
E+As, we find that M/L can be biased low by a factor 
of 1.25 to 5.5 relative to that derived using the red Tg. 

5.3.3. Will E+As Evolve into E/SOs? 

In summary, the scaling relations of E+As demon- 
strate that their stellar populations are currently dif- 
ferent from those of E/SOs. This conclusion is drawn 
from the FP parameters (cr, and /ig) in one photo- 
metric band, and therefore is derived independently from 
our spectroscopic knowledge that E-l-A galaxies have sig- 
nificant A-type populations. The surprising finding is 
that among E-|-As, the variation of stellar populations is 
closely tied to the structural parameters, i.e., E+As fol- 
low their own scaling relationships such that smaller or 
less massive galaxies have a smaller M/L. Such a trend 
arises naturally within a merger scenario, where low mass 
galaxies (the progenitors of low-mass E+As) have higher 
gas fractions (| Young fc Scovillj|1991h and could produce 
relatively larger populations of young stars. 

Will E+As fade onto the fundamental plane of E/SO 
galaxies after a few Gyr? This question is difficult to an- 
swer fully due to the complex interplay between galaxy's 
dynamical evolution and the evolution of its stellar pop- 
ulations. However, if we ignore the dynamical evolu- 
tion, we can estimate the amo unt of fading using a stel- 
lar population synthesis model (|Bruzual fc Charlotll2003l . 
BC03). We adopt a simple star formation history for 
the E+As: a single instantaneous starburst on top of 
an underlying 10 Gyr old single burst population. Both 
populations have Salpeter IMFs and solar metallicities. 

The largest uncertainties in determining the amount of 
fading and the time required for settling onto the E/SO 
FP are the unknown post-burst age (the time elapsed 
since the burst) and burst strength (the fraction of the 
stellar mass produced during the starburst). Because 
we are unable to break the degeneracy between post- 
burst age and burst-strength for this sample (see, how- 
ever, Yang et al. 2008, in prep.), we test whether the 
required fading time is reasonable for the galaxies whose 
post-burst ages and burst strengths are consistent with 
the FP offset (1.42 ± 0.36 mag in Gunn r). For vari- 
ous post-burst ages and burst fractions in Figure [T31 we 
calculate how much the galaxies should fade and how 
long it will take from the given post-burst age until their 
(B — R) colors agree with that of a 10 Gyr old simple stel- 
lar population, i.e., until they become E/SOs. The lined 
region in Figure [14] shows the post-burst ages and burst 
fractions that would satisfy the E+A selection criteria 
adopted by Zabludoff et al. (average Balmer equivalent 
width {H/3-fS) > 5. 5 A). The shaded region represents 
those galaxies that wiU fade by AM^ = 1.42 (±0.36) 
mag until their {B — R) colors match the 10 Gyr old 
population i.e., those E+As whose M/L's are consistent 
with the observed FP offsets. This shaded region illus- 
trates the degeneracy between the post-burst ages and 
burst strength: old/strong and young/weak bursts both 
explain the observed FP offsets. We also show the re- 
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Fig. 14. — Time required for E+As to fade until their {B — R) 
colors agree with that of 10 Gyr old simple stellar population as a 
function of burst fraction and post-burst age. From left to right, 
the contours represent fading times of 3, 5 and 7 Gyr, respectively. 
The lined region represents galaxies that would be classified as 
E+As according to the Zabludoff et al. (1996) selection criteria 
(average Balmer equivalent width (Hf3-yS) > 5.5A). The shaded 
region shows galaxies that will fade by their observed offsets from 
the E/SO fundamental plane {AMr = 1.42+0.36 mag) for the given 
burst fraction and post-burst age. This fading time lies between 3 
and 7 Gyr, indicating that E+As with reasonable post-burst ages 
(100 — 700 Myr) and burst fractions (7 — 50%) can evolve into 
E/SOs within 5 Gyr. 

quired fading time (3, 5, and 7 Gyr) with contours. The 
shaded region falls between 3 and 7 Gyr, indicating that 
E+A galaxies with reasonable post-burst ages (100 - 700 
Myr) and burst fractions (7 - 50%) will settle onto the 
E/SO FP after - 5 Gyr. 

6. YOUNG STAR CLUSTERS 

The properties of star cluster systems in early type 
galaxies, e.g., the color bimodality, provide important 
clues for how the host galaxies formed and how they 
assembled their stellar content over tim e (see the re- 
cent review by iBrodie fc Stradeil |2006[ ). If many or 
all E+A galaxies result from galaxy-galaxy close in- 
teractions or merg ers, as suggested here and in previ- 
ous work (.Zabludoffet al.lll99j|;lYMig et al.|[200l]200p ; 
IBlake et al.l [200l iTran et all [2iOOaT2004l: IGotol [2005h . 
and evolve into E/SOs, as we argue in this paper, the 
photometric properties of their cluster systems should 
be consistent, to within evolutionary corrections, with 
those of early types. In the WFPC2 sample, we discov- 
ered young star clusters around four E+As'' ()Yang et all 
|2004[ ). These clusters are much brighter than Galactic 
globular clusters and have blue colors consistent with 
ages estimated roughly from the E+A galaxy spectra 
(< Gyr). Our new ACS observations of the remaining 
LCRS E+As also reveal a number of point-like sources 
surrounding the galaxies. In this section, we first present 
the evidence that these are newly-formed star clusters, 
then investigate whether their colors and luminosities are 
consistent with being created concurrently with the burst 
that generated the young stars in the E+As, and finally 

■* EAOIA, 02, 03, 04 
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Fig. 15. — Distribution of the number of point sources in ran- 
dom beams sampling the ACS fields of the five E+As {EA07, 08, 
11, 15, 18) with a significant (> 90% confidence level) excess of 
surrounding sources. In each field, the beam size corresponds to a 
20 kpc radius at the E+A's redshift. In each panel, arrows repre- 
sent the number of cluster candidates within this radius centered 
on the E+A. The overlayed dashed lino is the Poisson distribution. 

test whether these cluster systems can evolve into the 
globular cluster systems of present-day E/SOs. 

6.1. Detection and Photometry of Cluster Candidates 

To detect cluster ca ndidates, we use SExtractor 
(|Bertin fc Arnoutill996D on both the original and model- 
subtracted residual R band images (Fig. [2|). After test- 
ing various selection criteria and visually inspecting the 
detected point sources, we require that cluster candidates 
have at least 4 adjacent pixels whose flux in each pixel is 
at least Susky larger than the local sky value, where (Tsky 
is the root mean square (rms) of the local sky values es- 
timated using a background mesh size of 4-8 pixels in 
SExtractor. We use the residual images to find compact 
sources within the galaxies, where the rapidly varying 
galaxy surface brightness prevents the reliable detection 
of sources in the original images. We choose to lower the 
detection threshold to 2.5 CTsky for the residual images, 
because dsky is inflated by the residual galaxy light. From 
this initial source list, objects with an apparent magni- 
tude niR < 22, elongation e > 1.7, and isophotal area 
^ > 36 pixels are rejected as possible foreground stars 
or background galaxies. 

Because clusters are expected to be unresolved at 
the distance of our E-l-A galaxies (~ 100 pc per pixel 
at z ^ 0.1), contamination from background galaxies, 
foreground stars, hot pixels, and residual cosmic-rays 
could be problematic. Therefore, we run a statistical 
test to determine whether there is an excess of com- 
pact sources near each E-|-A. First, we build a control 
sample of cluster-like sources from the entire ACS field 
(-200" X 200") outside of a 20 kpc radius from the F.+A 
using the same selection criteria. Then we calculate the 
distribution of the number of the sources within circular 
beams of 20 kpc projected radius at random positions in 
the ACS field outside a 20 kpc radius from the E-|-As. 
If there are significantly more cluster candidates within 
the 20 kpc beam centered on the E-t-A than in the ran- 
dom beams, then it is likely that these arc true compact 
objects associated with the E+A. We test various beam 
sizes, e.g., fixed angular size versus physical size, and find 



that our conclusions are insensitive to this choice. 

Five E+As (EA07, 08, 11, 15, 18) have an excess of 
associated point sources at a confidence level greater 
than 90%. We show the number of cluster candi- 
dates found near the E+As relative to the distribu- 
tion found within random beams in Figure 1151 Be- 
cause of the post-starburst nature of E+As and the 
merger signatures clearly seen in EA07, 11, and 18, it 
is likely that these sources are the newly formed star 
clus ters often found in on-going mergers (The Anten- 



nae; IWhitmore fc Schweizer 



fNaC7272: lWhitmore et~ 



19931 ) 



1995) and merger remnants 



To test whether the colors and luminosities of these 
cluster candidates are consistent with a merger origin 
and with evolution into the globular cluster populations 
of E/SOs, we measure aperture magnitudes in the B and 
R band images. Because many cluster candidates are 
located within the galaxy, where there is a rapidly vary- 
ing effective background, obtaining reliable photometry 
is difficult. We first subtract the background galaxy light 
by fitting 7th order polynomials to the background, ex- 
cluding the cluster itself, within a 11x11 pixel box cen- 
tered on the point source. We then add a fiat background 
that has the same mean as the fitted background to pre- 
serve the noise properties. Finally, we measure magni- 
tudes within a 2 pixel radius. Because the uncertainties 
in the measured magnitudes and colors transfer directly 
to those in the derived ages, we estimate photometric 
errors using artificial cluster tests. We place one thou- 
sand model star clusters with known brightness into the 
relevant parts of the images and repeat our procedure. 
We find that if the model cluster is detected, the mag- 
nitude error due to the background subtraction is nearly 
insensitive to the underlying average sky value and de- 
pends mostly on the brightness of the cluster. Typically, 
the magnitude uncertainty is ~0.25 mag down to ma — 
26.5. 



6.2. Luminosity Functions and Ages of the Star 
Clusters 

In Figure lTHl we present the luminosity functions (LFs) 
of cluster candidates in the four E+As (EAOlA, 07, 15, 
18) that have more than five clusters with measurable 
colors {(Jb-r < 0.35 mag) within a 30 kpc radius^ . The 
LFs of Milky Way globular clusters (iHarrid Il996l) and 
clu sters in a merging galaxy (NGC 3597: ICarlson et al.l 
1 1999) are also shown in the top panel for comparison. 
The thick solid lines and the shaded histograms repre- 
sent all E+A cluster candidates within 30 kpc and only 
those with measured colors, respectively. The detected 
clusters in the four E+As are generally much brighter 
than Galactic globular clusters (M/j < —11), but con- 
sistent with the bright end of the NGC3597 cluster LF, 
suggesting that these compact point sources are a popu- 
lation of star clusters formed during the starburst. 

More interesting is that the bright end of the LFs be- 
comes fainter as the EA number increases, i.e., as the 
host E+A becomes redder. Furthermore, we are able 
to detect a cluster population more often in E+As with 

^ Wc adopt a 20 kpc radius for the random beam test to deter- 
mine if cluster candidates exist, but we analyze all cluster candi- 
dates within a 30 kpc radius to include as many cluster candidates 
as possible. 
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Fig. 16. — Cluster luminosity functions for the four E+As 
(EAOl, 07, 15, 18) that have more than five clusters with measur- 
able colors. The bluest (and perhaps youngest) galaxy is at the 
top, the reddest galaxy at the bottom. The solid lines are the 
LPs for the cluster candidates projected within 30 kpc of the host 
galaxy. These LPs are incomplete fainter than about the middle of 
their range. The shaded histograms are for the clusters with mea- 
sured colors. The LPs of Galactic GCs and the merging galaxy 
NGC3597 are represented with dashed and dot-dashed lines, re- 
spectively. The clusters in the E+As are much brighter than the 
Galactic globular clusters (Mn < —11) but consistent with the 
bright end of the NGC3597 cluster LP, suggesting that these com- 
pact point sources are the population of star clusters formed during 
the starburst. The luminosity of the brightest clusters decreases as 
the color of the host galaxy reddens, suggesting a fading trend with 
age. The shaded regions are the globular cluster LPs of E/SOs that 
are evolved backward to the epoch of the E+A phase. The LPs 
of the cluster candidates fall in the shaded region, indicating that 
they could evolve into the globular cluster population of early-type 
galaxies. 

smaller EA number. We can reject with ^90% confi- 
dence that the distribution of E-|-As with cluster pop- 
ulations (EAOIA, 02, 03, 04, 07, 08, 11, 15, 18) is ran- 
dom. All nine of these E-f As exhibit tidal features or dis- 
turbed morphologies f ijA.ip . These trends suggest that 
the newly-formed young star cluster systems fade or dis- 
rupt as the merger remnant ages and morphologically re- 
laxes. If E-t-As evolve into E/SOs, we expect these cluster 
systems to become part of the globular cluster systems 
of E/SOs. 

To pursue these suggestions, we determine the ages of 
individual clusters by comparing the cluster colors with 
predictions from stellar population synthesis models. In 
Figure [T71 we show the {B — R) evolution for star clus- 
ters derived from a BC03 model and the colors of the star 
clusters found in three E -|-As (EA07 , 15, 18). We also 
reproduce Figure 10 from lYang et al.l (|2004D for the clus- 
ters in EAOIA for completeness. To avoid uncertainties 
in the K-correction and the magnitude transformation 
due to the unknown SED, we calculate the redshifted 
(S435— i?625) color evolution using the spectra provided 
by the BC03 stellar synthesis model. We assume an 
instantaneous stellar population with a Salpeter initial 
mass function (IMF) over stellar masses ranging from 1 
to 100 Mq and solar metallicity. On the model tracks 



in Figure [TTl we plot the observed cluster colors and the 
possible age ranges due to the photometric errors (shaded 
region) . We include only clusters with color uncertainties 
less than 0.35 mag. The ages of individual clusters are 
consistent with what one expects from the E-l-A spectra: 
older than ~10 Myr (the lack of emission lines associ- 
ated with massive OB stars) but younger than ^ Gyr 
(the presence of strong Balmer lines due to A stars). 

To use the clusters to derive the time since the star- 
burst, when the majority of the clusters and the young 
galactic stars formed, we apply a statistical test assum- 
ing that all star clusters formed during a single burst 
and therefore that their spread in observed color arises 
only from measurement errors. Details of the post-burst 
age estimation and statistical tests of the coeval assump- 
tion are described in Appendix [BJ We show the post- 
burst age inferred from the clusters as horizontal bars in 
Figure [171 EAOIA, 07, 15, and 18 have post-burst age 
ranges of [10 Myr, 450 Myr], [100 Myr, 1.5 Gyr], [150 Myr, 
IGyr], and [400 Myr, IGyr], respectively, corresponding 
to the 95% confidence level. Note that it may be pos- 
sible to break the long-standing degeneracy between the 
burst strength and the time elapsed since the starburst 
by using individual cluster ages, because these systems 
are much simpler than the composite population of the 
host galaxies. However, small number statistics and low 
S/N in the color measurements do not allow us to put 
tighter constraints on the post-burst ages at this time. 

The evolution of the cluster LFs is consistent with the 
model predictions. While the uncertainties are fairly 
large, the derived post-burst ages are consistent with the 
E-t-A spectra and tend to increase with EA number, sug- 
gesting that the EA numbers are roughly consistent with 
a sequence of post-burst ages. In Figure[T8l we show that 
the evolution of the bright magnitude limit of the cluster 
LF and the post-burst ages (derived from the cluster col- 
ors) are consistent with the fading of a uniform cluster 
mass function. 



6.3. Evolution of the Cluster Population 

The remaining question is whether these young cluster 
candidates can evolve into a significant part of the glob- 
ular cluster systems of present-day E/SOs. Specifically, 
are the current numbers of clusters in E-|-As consistent 
with the globular cluster luminosity function (GCLFs) 
of E/SOs? It has been proposed that the metal- rich 
GCs in giant ellipticals formed during mergers of gas-rich 
spirals, producing the bimodality of the GC subpopula- 
tions (|Ashman fc Zepilll992l) . In deed, many youiig GC s 
hav e been found in mergers (e.g.. lWhitmore et al.lll993[) . 
and lGoudfrooii et al.l (|2004D demonstrated that the LF of 
intermediate-age (~3 Gyr) GCs in NGC1316 can evolve 
dynamically into the red metal-rich cluster populations 
that are ubiquitous in E/SOs. Unfortunately, due to low 
S/N in the color measurements and the small number 
statistics of our cluster candidates, it is difficult to pre- 
dict the evolution of the individual cluster systems in our 
sample. Instead, we "predict" the cluster LFs of E-|-As 
under the assumption that they are indeed the progeni- 
tors of the E/SO cluster population. 

First, we assume the following form of the GCLFs, i.e.. 
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Fig. 17. — Evolution of the redshifted (B — R) colors of star 
clusters derived from BC03 and th e clu ster colors from the four 
E+As (EAOl, 07, 15, 18) in Figure [16] For the model track, we 
assume an instantaneous stellar population with a Salpeter IMF, 
masses ranging from 1 to 100 Mq, and solar metallicity. On the 
model track, wo show the measured cluster colors and the possible 
range of ages due to the photometric errors (shaded region). Only 
the clusters with color uncertainties less than 0.35 mag are included 
in the age determination. The ages of clusters derived from their 
colors are consistent with our expectations from the E+A galaxy 
spectra: older than ~10 Myr (the lack of the emission lines due to 
the massive OB stars), but younger than ~ Gyr (the presence of 
strong Balmer lines due to the A stars). The arrows represent the 
estimated post-burst age (time elapsed since the starburst) using 
the methods described in i|6.2l Cluster ages become older as the 
EA number increases, i.e., as the host E+A becomes redder. 
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where My is the total F-band magnitude of the host 
galaxy, Mq and a are the turnover magnitude and 
the dispersion of the Gaussian function, Sn is the 
specific frequency given by S n = Nqq ■ iQ^-'^i'^'^v+is) 
(jHarris &: van den Berghlll981), and / mr is the fraction 
of metal-rich "red" clusters. iHarrisI (|2001[ ) found that 
(Sn) = 3.6 ± 0.5, Mq = -7.4 ± 0.15 mag in the F-band 



and a = 1.4 ± 0.03 for local elliptical galaxies {My < 
-18). 

Second, to evolve this GCLF backward to the E-l-A 
epoch, we calculate the fading of the turnover magnitude 
between an age corresponding to the E-f A clusters and 
that corresponding to E/SG GCs, assuming that E/SO 
GCs are 5-10 Gyr old. For solar metallicity and a 
Salpeter IMF, a BC03 model predicts ~ 2 mag of fading 
between 1 Gyr and 10 Gyr. Note that these derived LFs 
will be lower limits because we consider only photometric 
evolution. For host galaxy luminosities, we adopt My 
— My^'^ + (1.3 ± 0.73), as measured from the scaling 



-15 



-14 



T3 
3 



-13 



-12 



-11 



v5xl0°tf, 



\ 



. 2xlO°Af„ 



vl0''Mo\_. 



lJ_ 



10 



100 

Burst Age [Myr] 



1000 



Fig. 18. — Evolution of the bright end of t he st ar clus ter LFs in 
the four E+As (EAOIA, 07, 15, 18) in Figures [16] and [iTl We show 
the relation between the average luminosity of the five brightest 
clusters in each E-l-A and the post-burst age shown in Figure [T71 
Dot-dashed lines represent BC03 model tracks for simple stellar 
populations (SSPs) with different cluster masses (1, 2, 5 and 10 X 
10^ Mq). The observed fading of the cluster LFs is consistent with 
the color-derived post-burst ages. 

relationships. 

The "predicted" GCLFs taking into account all possi- 
ble errors mentioned above are consistent with the ob- 
served E-l-A cluster LFs (Figure [TB)) . The most uncer- 
tain parameter is the age difference between the clusters 
in E-|-As and in E/SOs due to the fast luminosity evo- 
lution around 1 Gyr. We use ages of 100 - 500 Myr 
and 500 - 1000 Myr for EAOlA and the other E-f As, 
respectively. We vary /mr between i and |. Most of 
our observed LFs fall within the uncertainties (shaded 
regions), indicating that the number of clusters observed 
does not violate the constraint derived from the GCLFs 
of the present-day E/SOs. While it is highly uncertain 
how many clusters in E-|-As will survive and how much 
they will contribute to the metal-rich GCs in E/SOs, we 
conclude that it is at least possible that young star clus- 
ters in E-|-As can evolve into the globular cluster systems 
of E/SOs. 

7. CONCLUSIONS 

We study the detailed morphologies of 21 E-l-A 
galaxies using high resolution HST/ACS and WFPC2 
images to investigate into what E-l-A galaxies will evolve 
after their young stellar populations fade away in a few 
Gyr. Our findings are: 

1. The morphologies of E-|-As are extremely diverse, 
ranging across train-wrecks, barred galaxies, and blue- 
cores to relaxed disky galaxies. Most of these galaxies 
lie in the field, well outside of rich clusters, and at least 
11 (55%) have tidal or other disturbed features. Our 
sample includes one binary E-l-A system, in which both 
E-fAs are tidally disturbed and interacting with each 
other. These results support the picture in which galaxy- 
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galaxy tidal interactions or mergers are responsible for 
triggering the E+A phase in many cases. 

2. E+As are bulge-dominated systems (median bulge 
fraction B/T = 0.59) and their light is highly concen- 
trated (Sersic index n > 5). When dust is negligible (at 
least 67% of the time), E-|-As have high concentration 
indices (C > 4.3) consistent with those of spheroids, but 
considerably larger asymmetry indices {A > 0.04) than 
ellipticals due to structures within a few that presum- 
ably arise from the starburst and/or recent merger. Thus 
E-|-As would be morphologically classified as early-type 
galaxies once these disturbances relax and the low sur- 
face brightness tidal features dissipate or fade. 

3. The color morphologies of E-|-As are as diverse as 
their structural morphologies. A large fraction (70%) 
have positive color gradients (bluer toward center), in- 
dicating that their young stellar populations are more 
concentrated than their older populations. We demon- 
strate that evolution can invert these gradients into the 
negative gradients typical of E/SOs if the inner parts of 
E-|-As have become more metal enriched than the outer 
parts due to the centralized star formation. 

4. Six E+As (30%) exhibit compact (0.4-1.4 kpc) blue 
cores, which mig ht be the local analogs o f the high-z ellip- 
tical blue-cores ()Menanteau et al.ll2001[ ). W e discovered 
LINE Rs in three of these blue-core E-|-As (I Yang et al.l 
l2006l ). and the relationship between LINERs and blue- 
cores could be an important clue to what stops the star 
formation in E-l-A galaxies. 

5. E-|-As stand apart from the E/SO fundamental plane 
(FP) in the edge-on projection, implying that the stellar 
populations of E-|-As are different from that of E/SOs. 
E-|-As have, on average, a M/L that is 3.8 times smaller 
than that of E/SOs. The tilt of the E+A FP indicates 
that the variation of the stellar populations among E-|-As 
is closely tied to their structural parameters, i.e., E-t-As 
follow their own scaling relationships such that smaller 
or less massive galaxies have smaller M/L. Such a trend 
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arises naturally within a merger scenario, where low mass 
galaxies (the progenitors of low-mass E-|-As) have higher 
gas fractions (| Young fc Scovilldl 19911 ) and could produce 
relatively larger populations of young stars. 

6. We find a population of unresolved compact sources 
in at least nine E-f A galaxies (45%). The colors and 
luminosities of these young star cluster candidates are 
consistent with the ages inferred from the E-l-A spectra 
(0.01 - 1 Gyr). The bright end of the cluster luminosity 
function fades as the host galaxy becomes redder, sug- 
gesting that the newly-formed young star clusters age in 
parallel to their host. This interpretation is confirmed 
by the color evolution of the cluster systems. 

We have now examined the full set of E-|-As from the 
Las Campanas Redshift Survey and so have representa- 
tive results for local E-f A galaxies. We have used high 
spatial resolution images to probe their detailed mor- 
phologies. We find that their properties are either con- 
sistent with those of E/SOs or, if left to evolve passively, 
will become like those of early-types. The morphologies, 
color profiles, scaling relations, and young star clusters 
suggest that E-|-As galaxies are caught in the act of trans- 
forming from late- type to early- type galaxies. 



We thank an anonymous referee for his/her thorough 
reading of the manuscript and helpful comments. We 
thank Chien Peng for help in using the GALFIT and use- 
ful discussions about galaxy morphology. YY and AIZ 
acknowledge support from NSF grant AST-02-06084, 
and ffST grant GO-9781. Support for proposal GO-9781 
was provided by NASA through a grant from the Space 
Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, 
Incorporated, under NASA contract NAS5-26555. 

Facilities: HST(ACS,WFPC2) 



Chang, T., van Gorkom, J. H., Zabludoff, A. I., Zaritsky, D., & 

Mihos, J. C. 2001, AJ, 121, 1965 
Conselice, C. J., Bershady, M. A., & Jangren, A. 2000, ApJ, 529, 

886 

Couch, W. J. & Sharpies, R. M. 1987, MNRAS, 229, 423 

Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 

Dressier, A. 1987, ApJ, 317, 1 

Dressier, A. & Gunn, J. E. 1983, ApJ, 270, 7 

Faber, S. M., & Jackson, R. E. 1976, ApJ, 204, 668 

Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 

Franx, M. 1993, ApJ, 407, L5 

Franx, M., lUingworth, G., & Heckman, T. 1989, AJ, 98, 538 
Frei, Z., Guhathakurta, P., Gunn, J. E., & Tyson, J. A. 1996, AJ, 
111, 174 

Fukugita, M., Shimasaku, K., & Oshawa, T. 1995, PASP, 107, 945 

Galaz, G. 2000, AJ, 119, 2118 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Goto, T. 2004, A&A, 427, 125 

Goto, T. 2005, MNRAS, 357, 937 

Goto, T., et al. 2003, PAS J, 55, 771 

Goudfrooij, P., Gilmore, D., Whitmore, B. C, & Schweizer, F. 

2004, ApJ, 613, L121 
Graham, A. W. 2002, MNRAS, 334, 859 

Graham, A. W., Erwin, P., Caon, N., & Trujillo, I. 2001, ApJ, 563, 
Lll 

Gunn, J. E., & Gott, J. R. I. 1972, ApJ, 176, 1 

Hamabe, M., & Kormendy, J. 1987, lAU Symp. 127: Structure and 

Dynamics of Elliptical Galaxies, 127, 379 
Harris, W. E. 1996, AJ, 112, 1487 



26 



Yang et al. 



Harris, W. E. 2001, Saas-Fee Advanced Course 28: Star Clusters, 
223 

Harris, W. E., & van den Bergh, S. 1981, AJ, 86, 1627 

Hogg, D. W., Masjedi, M., Berlind, A. A., Blanton, M. R., 

Quintero, A. D., & Brinkmann, J. 2006, ApJ, 650, 763 
Jedrzejewski, R. I. 1987, MNRAS, 226, 747 

j0rgensen, I., Franx, M., & Kjaergaard, P. 1996, MNRAS, 280, 167 
Kelson, D. D., lUingworth, G. D., van Dokkum, P. C, & Franx, 

M. 2000, ApJ, 531, 137 
Kelson, D. D., lUingworth, G. D., van Dokkum, P. C, & Franx, 

M. 2000, ApJ, 531, 184 
Kennicutt, R. C. 1998, ARA&A, 36, 189 
Kent, S. M. 1985, ApJS, 59, 115 
Kobayashi, C. 2004, MNRAS, 347, 740 
Kormendy, J. 1977, ApJ, 218, 333 

Lee, J. H., Lee, M. G., & Hwang, H. S. 2006, ApJ, 650, 148 
Lotz, J. M., Primack, J., & Madau, P. 2004, AJ, 128, 163 
Lupton, R., Blanton, M. R., Fekete, C, Hogg, D. W., O'MuUane, 

W., Szalay, A., & Wherry, N. 2004, PASP, 116, 133 
Menanteau, F., Abraham, R. C, & EUis, R. S. 2001, MNRAS, 322, 

1 

Menanteau, F., et al. 2004, ApJ, 612, 202 
Mihos, J. C, & Hernquist, L. 1994, ApJ, 425, L13 
Mihos, J. C, & Hernquist, L. 1994, ApJ, 427, 112 
Miller, N. A. & Owen, F. N. 2001, ApJ, 554, L25 
Norton, S. A., Gebhardt, K., Zabludoff, A. L, & Zaritsky, D. 2001, 
ApJ, 557, 150 

Pahre, M. A., Djorgovski, S. G., & de Carvalho, R. R. 1998, AJ, 
116, 1591 

Peletier, R. F., Davies, R. L., lUingworth, G. D., Davis, L. E., & 

Cawson, M. 1990, AJ, 100, 1091 
Peng, C. Y., Ho, L. C, Impey, C. D., & Rix, H. 2002, AJ, 124, 266 
Petrosian, V. 1976, ApJ, 209, LI 

Poggianti, B. M., Small, I., Dressier, A., Couch, W. J., Barger, 
A. J., Butcher, H., Ellis, R. S., & Oemler, A. J. 1999, ApJ, 518, 
576 

Pracy, M. B., Couch, W. J., Blake, C, Bekki, K., Harrison, C, 
CoUess, M., Kuntschner, H., & de Propris, R. 2005, MNRAS, 
359, 1421 

Quintero, A. D., et al. 2004, ApJ, 602, 190 
Rix, H. & Zaritsky, D. 1995, ApJ, 447, 82 

Rose, J. A., Gaba, A. E., CaldweU, N., & Chaboyer, B. 2001, AJ, 
121, 793 



Rothberg, B., & Joseph, R. D. 2006, AJ, 131, 185 

Schade, D., Lilly, S. J., Crampton, D., Hammer, F., Le Fevre, O., 

& Tresse, L. 1995, ApJ, 451, LI 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Shectman, S. A., Landy, S. D., Oemler, A., Tucker, D. L., Lin, H., 

Kirshner, R. P., & Schechter, P. L. 1996, ApJ, 470, 172 
Sirianni, M., et al. 2005, PASP, 117, 1049 

Small, L, Morrison, G., Gray, M. E., Owen, F. N., Ivison, R. J., 

Kneib, J.-P., & Ellis, R. S. 1999, ApJ, 525, 609 
Tacconi, L. J., Genzel, R., Lutz, D., Rigopoulou, D., Baker, A. J., 

Iserlohe, C, & Tecza, M. 2002, ApJ, 580, 73 
Tamura, N., Kobayashi, C, Arimoto, N., Kodama, T., & Ohta, K. 

2000, AJ, 119, 2134 
Tran, K.-V. H., Franx, M., lUingworth, G., Kelson, D. D., & van 

Dokkum, P. 2003, ApJ, 599, 865 
Tran, K.-V. H., Franx, M., lUingworth, G. D., van Dokkum, P., 

Kelson, D. D., & Magee, D. 2004, ApJ, 609, 683 
Treu, T., et al. 2005, ApJ, 633, 174 

Treu, T., Koopmans, L. V., Bolton, A. S., Buries, S., & Moustakas, 

L. A. 2006, ApJ, 640, 662 
van Dokkum, P. G. 2001, PASP, 113, 1420 
van Dokkum, P. G., & EUis, R. S. 2003, ApJ, 592, L53 
Whitmore, B. C, Schweizer, F., Leitherer, C, Borne, K., & Robert, 

C. 1993, AJ, 106, 1354 
Whitmore, B. C, & Schweizer, F. 1995, AJ, 109, 960 
Wise, M. W., & Silva, D. R. 1996, ApJ, 461, 155 
Yagi, M., Goto, T., & Hattori, T. 2006, ApJ, 642, 152 
Yamauchi, C, & Goto, T. 2005, MNRAS, 359, 1557 
Yan, R., et al. 2008, ArXiv e-prints, 805, arXiv:0805.0004 
Yang, Y., Zabludoff, A. L, Zaritsky, D., Lauer, T. R., & Mihos, 

J. C. 2004, ApJ, 607, 258 
Yang, Y., Tremonti, C. A., Zabludoff, A. I., & Zaritsky, D. 2006, 

ApJ, 646, L33 

Yang, Y., Tremonti, C. A., Zabludoff, A. I., 2008, in preparation. 
Young, J. S., & Scoville, N. Z. 1991, ARA&A, 29, 581 
Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 803 
Zabludoff, A. L, Zaritsky, D., Lin, H., Tucker, D., Hashimoto, Y., 

Shectman, S. A., Oemler, A., & Kirshner, R. P. 1996, ApJ, 466, 

104 

Zaritsky, D., Gonzalez, A. H., & Zabludoff, A. I. 2006, ApJ, 638,725 



TABLE 1 
Properties of E+A Galaxies 



ID'' 




R.A. 


Dec. 


2=^ 


Scale 


^ b 


SFR 


H I mass 


Morphology 




(J2000) 


(J2000) 




(kpc/") 


(km s-l) 


(Me yr-i) 


(109 Mq) 




EAOIA 


11 


01 19.05 


1 o in 1 '7 o 


U.U/4DO 


1.42 


Qc+39 
■^■^-35 


< (J. 12 ± U.U^ 


T '7 J- n /if 
/.7 ± U.4 


tidal-feature, companion, cluster 


EAOIB 


1 1 
i 1 


Ui io.oy 


— 12 10 14.1 




1.42 


^ n n7 _|_ n hqc 

^ yj.KJl _l_ yj.KJrj 




1 1 H n 1— Tpn 1 11 rfi rTimr^amr^n T.TlMTi.T-? 

LlLldl 1 Kdli Lli t; , ^jUlllUtXillUli, 1 j 1 1 \ 1 1 1 1 


EA02 


02 


1 7 as 


-44 32 47.2 


0.09873 


1.82 


202+}! 
^ — 16 


< 4.4 


< 3.4^ 


tidal-feature, dust?, cluster 


EA03 


12 


09 5.54 


-12 22 35.8 


0.08109 


1.53 


120+on 
— 20 


< 1.7 " 


< 4.2^ 


tidal-feature, cluster 


EA04 


04 


00 0.34 


-44 35 15.8 


0.10122 


1.86 


131+^ 

^ — 9 


< 3.2 


< 2.3^ 


disturbed, dust?, cluster 


EA05 


01 


58 1.52 


-44 37 14.0 


0.11731 


2.12 


120^8 






bar, blue-core 


EA06 


11 


^'^ 


-03 10 36.3 


0.08850 


1.65 


^^ — 23 


< 0.10 ± 0.03'' 




dust, blue-core, LINER 


EA07 


22 


41 9.75 


-38 34 34.9 


0.11413 


2.07 


241+U 

— 12 






tidal-feature, dust, cluster 


EA08 


14 




-12 57 33.1 


0.11222 


2.04 


99+?, 


< 1.8 ° 




tidal-feature, cluster 


EA09 


01 


J- / OO . ^o 


-41 34 24.6 


0.06513 


1.25 


74+134 
' —5 


< 1.3 " 




blue-core 


EAIO 


02 


11 42.96 


-44 07 40.2 


0.10493 


1.92 


121+8 






tidal-feature, companion 


EAll 


01 


14 49.65 


-41 22 31.8 


0.12167 


2.19 


174+;? 

' —11 


< 9.6 


< 5.3^ 


disturbed, companion, cluster 


EA12 


12 


05 59.79 


-02 54 32.4 


0.09713 


1.80 


111+;^ 


7.9 <= 




tidal-feature, companion, dust 


EA13 


11 


19 52.44 


-12 52 39.6 


0.09572 


1.77 


165+10 


< 2.1 




tidal-feature, companion, dust 


EA14 


13 


57 1.68 


-12 26 48.2 


0.07046 


1.34 


135li° 


< 1.5 " 




bar, blue-core 


EA15 


14 


40 44.26 


-06 39 52.3 


0.11381 


2.07 


125l?o 


< 0.05 ±0.02'^ 




cluster 


EA16 


12 


19 55.69 


-06 14 1.0 


0.07642 


1.45 




< 0.05 ±0.02'^ 




blue-core, LINER 


EA17 


10 


13 52.40 


-02 55 47.2 


0.06090 


1.18 


23+1 


< 0.06 ± 0.03'^ 


< 2.9s 


dust, blue-core, LINER 


EA18 


00 


22 46.79 


-41 33 35.8 


0.05985 


1.16 


23+" 


< 1.2 "= 


2.0 ± 0.3s 


disturbed, dust, cluster 


EA19 


02 


07 49.45 


-45 20 51.1 


0.06400 


1.23 


84+18 
"^-34 


2.7 " 


< 1.2s 


dust 


EA21 


11 


15 24.96 


-06 45 13.8 


0.09944 


1.84 


1341^ 


< 1.7 <= 




bar 



E-l-A IDs and redshifts from Zabludoff ct al. ( 1996).'' Velocity dispersion of the old stellar population from lNorton et al.i (|2001i l. The entries with ct = 23 km s 

) = Lhq, / 1.27 X IQ-iygs s ~i (iKennicuttlliggl yi SFR from 
(W Hz"i) (|Yun et al.ll2001l '). Corrected for the cosmological 
' ' SFR from the 



are the upper limits. SFR from Ha luminosity ( Yang et al.|[200gl using the relation, SFR (M( 



yr 



yr 



= 5.9 X 10" 



the radio-continuum (jChang et al.ll200ll ) using the relation, SFR {M(_ 

parameters adopted in this paper. Note that we do not adopt the SFRs de rived from th e gas surface densities, which are highly uncertain 
radio-continuum (Miller & Owen 20 0l|). Same as footnote d.* H I mass from [Chang et ahl ()2001l '). Corrected for the cosmological parameters adopted in this 
paper. s H I mass from iBuvle et al.i ((20061 ') . 



+ 
> 

E» 
P 

H 

Q 

X 



to 
-J 



Yang et al. 



TABLE 2 
Magnitude Transformation 
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.63 


± 


0, 


,06 


3, 


,84 


± 


0, 


,15 


0, 


,04 


18, 


,47 


± 


0, 


,02 


17, 


,30 


± 


0.04 


EA17 


2 


.84 


± 


0, 


.05 


5. 


,04 


± 


0, 


.01 


0, 


,08 


18, 


,58 


± 


0, 


.01 


17, 


,16 


d= 


0.01 


EA18 


2 


.92 




0, 


.05 


3. 


,73 


± 


0, 


.01 


0, 


,19 


17, 


,81 


± 


0, 


.01 


16, 


,37 


± 


0.01 


EA19 


4, 


.43 




0, 


.03 


2. 


,95 




0, 


.01 


0, 


.10 


18, 


,06 


± 


0, 


.01 


16, 


,66 


± 


0.01 


EA21 


1, 


.22 




0, 


.17 


4. 


,64 


± 


0, 


.18 


0, 


.04 


18, 


,85 


± 


0, 


.04 


17, 


,26 


± 


0.05 



Typical errors for A are a{A)/A ~ 0.015. 
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TABLE 4 

DE VAUCOULEURS AND SeRSIC PROFILE FiT {R702 OR ^626 BANDS) 



















Scrsic 


T'-'^/^-Law 










ID 










c 


xi 




Me 


n 






c 


xi 




(kpc) 


(mag arcsec"-^) 


(fe/a) 






(kpc) 


(mag arcsec~2) 




(b/a) 






EAOIB 




08+°- 


^^gg+0.12 


0.70 


—0.04 


2.88 


1 35+°-*° 
^•■'"^ -0.01 


18 ,39 +0-55 


5 3g+0.99 
"^■"^^ -0.05 


0. 


.68 


—0.06 


2.53 


EA02 


2 


■'"-0.01 


19 83+°" 


0.73 


0.15 


1.51 


5 85+^°' 

-0.01 


21 41+°" 
^^•^^ -0.35 


6 79+°-™ 

-0.48 


0. 


.75 


0.34 


1.43 


EA03 


I 


fiq+0-09 
-0.01 


18 32+°-^° 
-0.05 


0.97 


—0.04 


2.26 


^■■^■^ -0.01 


20 65 +0-3'i' 
.iu.uu _o 22 


9 55+°-'^'' 

"■■J^ -0.37 


Q 


98 


—0.04 


1.50 


EA04 


2 


19+0.11 
•^^-0.01 


18 92+°-^° 


0.97 


0.01 


1.42 


-0.01 


19 76 +0-23 
-0.11 


5 37+0.29 

-0.17 


0. 


.98 


—0.02 


1.31 


EA05 


1, 


.62+°- 

" —0.01 


18 57+° °'' 

10.01 _Q 05 


0.61 


—0.10 


1.37 


9 Qo +0.82 
^•^•^ -0.01 


19 91+0-51 


7Qg+0.89 
'■^"-0.41 


0. 


.62 


—0.06 


1.30 


EA06 


Q 


40 +0.01 
■^■^ -0.01 


17 07+° °^ 

'^'■"1 -0.04 


0.86 


—0.05 


0.70 


+0.01 
-0.01 


17 25+0-05 

-^'■^■^-0.07 


4 60+°-'^ 


0. 


.86 


—0.05 


0.69 


EA07 


11 


7. +0.19 

•'°-o.oi 


21 27+°°'' 

-0.05 


0.74 


0.07 


1.27 


-07+0.35 

-0.01 


20 44+°-"* 


2 72+0.11 
^- ' ^ -0.09 


0. 


.75 


0.07 


1.12 


EA08 


1 


84+0-05 
■"*-o.oi 


19.58 +°°f 

—0.02 


0.89 


0.07 


0.99 


4.00+°-^? 

—0.01 


21.25+°?^ 
—0.19 


7.67+;;-^;' 


0. 


.89 


0.08 


0.85 


EA09 


2 


10 +0.08 
■'-^ -0.01 


20 86+° °® 


0.52 


—0.37 


1.61 


20 92 +-^2.12 


25 58+°-*'' 


15.22 


0, 


52 


—0.37 


1.39 


EAIO 


2 


99 +0.03 
•^^ -0.01 


19 73+0-02 
ia./o_o 05 


0.81 


0.05 


0.83 


3 28+0.70 
■^■^^ -0.01 


20 57+"-*' 


5 74+1^00 

°- -0.31 


0, 


,81 


0.04 


0.75 


EAll 


1, 


on +0.04 
•80-0.01 


18 92+° °* 


0.53 


0.14 


1.16 


2 8g+0.44 
-0.01 


19 97+0-29 

— t). 15 


6 31 

"■■^^ -0.24 


0. 


.53 


0.13 


1.06 


EA12 


4, 


g7+0.25 
-0.01 


21 47+°-°^ 


0.81 


1.13 


1.29 


g 14+2.40 
-0.01 


22 83+°-*' 

^^.00 _o 33 


Q 33+0.56 
" -0.40 


0. 


.81 


1.25 


1.22 


EA13 


5, 


-9+0.28 
•'^-0.01 


20.501°;- 


0.55 


0.22 


0.94 


8 52 +^•"1 
^•^^-0.01 


21 35 +°-'''^ 


5 56+°-'^2 
"'■^"-0.48 


0, 


,54 


0.17 


0.90 


EA14 


1 


^0+0.01 
■^^ -0.01 


,0 90+0.00 


0.57 


0.18 


1.68 


144+0-01 
-0.01 


18.83 +°,°i 


5.30 1°;- 


0, 


,56 


0.20 


1.58 


EA15 


1, 


45+0.01 

■^"^-0.01 


19.231°;°° 


0.91 


-0.07 


0.83 


9 00 +0.09 
^••^«-0.01 


20 34+o:o8 


6.54 1°;- 


0. 


,91 


-0.07 


0.75 


EA16 





+0.01 

•°^-o.oi 


18.041°;°^ 


0.83 


0.05 


1.25 


171+0.13 
^•'^-0.01 


20-18^0:31 


Q QO ~!"0 -33 
«-^^-0.70 


0, 


,81 


0.07 


0.84 


EA17 


1, 


„.+0.07 
■"^^ -0.01 


20 02 +0-08 


0.65 


-0.13 


1.28 


+3.89 

^■-i°-o.oi 


22.82+^,;,^" 


10 44 


0, 


,64 


-0.11 


0.96 


EA18 


3 


,-+0.03 
-0.01 


20 57+° °l 


0.50 


-0.27 


1.66 


2 24+0-04 
^•^"i-O.Ol 


19 76+0-04 
la. (D _o og 


2gg+0.06 

Z.DU _g og 


0, 


,49 


-0.29 


1.42 


EA19 


4, 


no +0.18 

■^°-o.oi 


21.411°:°^ 


0.43 


-0.14 


1.65 


Q 99+0.04 
•^•^^-0.01 


20 16+°°3 

ZU.ID _o 02 


lQg+0.03 
J^-"" -0.02 


0. 


,40 


-0.10 


1.06 


EA21 





Q9+0.01 

■^"'-n.oT 


18.20 1°- 


0.72 


0.21 


1.28 


2.20 l^n? 


20.091°-^® 


g op. +0.60 


0. 


,70 


0.19 


1.04 



Note. — EAOIA is excluded because it is so irregular. 
Typical errors for q and c are ~ 0.01. 



TABLE 5 

Bulge-Disk Decomposition (il702 and i?625 bands) 



rV*-Bulge Disk 



ID 


(kpc) 


Me 

(mag arcsec~2) 


(b/a) 


c 


rd 
(kpc) 


(mag arcsec~2) 


9^ 
{b/a) 






B/T 


xl 


EAOIB 


„ 73+0.01 

^J- '•5-0.01 


17 05+0.01 
^'•■5° -0.04 


0.75 


-0.06 


n QO ~l~0. 23 

^••^8-0.01 


19.351°;- 


0.39 


1, 


.21 


n 7S+0-01 

0.'8_o.03 


2.40 


EA03 


Q 71 +0.01 

' ^ -0.01 


16 94+0-03 
iD.ai _o 04 


0.93 


0.09 


Q cn +0.17 

3.60^0 01 


19.191°;- 


0.88 


0, 


,04 


54+0 01 
"•'-'^ -0.02 


1.57 


EA04 


1 11+0.01 
^■^^ -0.01 


1703+0.02 


0.87 


0.08 


3 40+0.11 
■J-^" -0.01 


19 30+° °*' 
la.ou 


0.72 


0, 


,32 


o.6il°:°3 


1.24 


EA05 


, ,9+0.04 
-^■^^ -0.01 


17 96+0-05 

-'-'•-^"-0.07 


0.57 


-0.15 


4 3P;+0.98 


20 59 +0--''2 


0.77 


-0, 


.38 


0.681°;- 


1.30 


EA06 


0.401°;- 


16 94 +"■"■' 

^ -0.03 


0.85 


-0.05 


, 09 +4.05 
4-8^_o.oi 


22-49 1°;S 


0.55 


0, 


,84 


n at; +0-04 
0.83 _„.()<) 


0.67 


EA07'' 


11.59 1°;« 


21 43+0-03 

^-^■^•5-0.07 


0.69 


0.17 


oc:+0-G3 
-0.01 


19.251°;- 


0.82 


-0, 


,50 


0.85 1";- 


1.09 


EA08 


0.59 1°;- 


1703+0.08 
i/.ao_o 00 


0.89 


-0.03 


2 16+°" 

^■^O-O.Ol 


19 24+° '° 


0.90 


0, 


.25 


0.45 1;;;- 


0.82 


EA09 


06+°°' 

-0.01 


15-641°;°: 


0.98 


-0.81 


145+0.02 
-0.01 


10 7-1 +0.02 
J^=-'^-0.00 


0.41 


-0, 


.21 


16+°°' 
"-0.01 


0.84 


EAIO 


1 16+°- 

-0.01 


18.901°;- 


0.85 


0.17 


,0+0.18 
•5-J^'5_o.oi 


1Q 9s +015 


0.68 


-0, 


,12 


64+°°' 
"."^-0.12 


0.74 


EAll 


lQg+0.03 
-0.01 


18.00 1°;- 


0.42 


0.47 


3 24+0.36 
•^■^^-0.01 


20.101°;- 


0.95 


-0, 


,68 


61+°°' 
^■"^ -0.02 


0.92 


EA12'= 


1 49 +°-''' 

-0.01 


19 74+0.10 

^«-''*-0.05 


0.65 


-0.48 


00 +0.56 
3.88 _o 01 


19.82 1°;« 


0.69 


-0, 


.56 


36 +°°"' 
U..3D_o oj 


1.04 


EA13 


-0.01 


19 24+° '® 

J^y-^''-0.25 


0.61 


0.10 


5 07+° ®' 
'^•"' -0.01 


19.431-1 


0.43 


0, 


,07 


47+0.05 

".^' -0.05 


0.83 


EA14 


75+0-01 
^■'^-0.01 


17.621°;°° 


0.55 


0.10 


9 00 +0.05 
'^■■^■^ -0.01 


19.881°;- 


0.40 


0, 


,33 


7+0.01 

" ' -0.01 


1.36 


EA15 


Q 47+0.01 

"•^' -0.01 


17.681°;- 


0.92 


0.02 


, 57+0.01 


18.761°;- 


0.90 


-0, 


,11 


n /IS +0-01 
0.48_o.o2 


0.73 


EA16 


Q 34+0.01 

"■■5^-0.01 


17.031°;- 


0.81 


0.02 


9 9, +0.09 
^•^^ -0.01 


20.02 1°;°» 


0.81 


0, 


,14 


57+0-01 
-0.01 


0.72 


EA17 


25+° °' 


17.481°;- 


0.69 


-0.15 


1 62+°°^ 
^■0^-0.01 


19.36+°°^ 


0.58 


0, 


,06 


35+0-01 

'^..5^-0.03 


0.78 


EA18'= 


9 Qo+0.02 
^•^'^-O.Ol 


20 59+°-°° 
zu.oa 0.05 


0.52 


-0.30 


63+°°' 
^•^"^-0.01 


18.641°:- 


1.00 


-0, 


,95 


91+°°' 
".^1 -0.01 


1.24 


EA21 


„ +0.01 


17 41 +° °^ 

^'•^^ -0.03 


0.70 


0.33 


4ig+0.32 
^•^^ -0.01 


20.60 1°;1° 


0.65 


-0, 


,23 


59+°- 

".3^-0.02 


0.93 



Note. — EAOIA, EA02 and EA19 arc excluded from the bulge-disk decomposition because the GALFIT program 
fit does not converge. 

^ Typical errors for q and c are ~ 0.01.'' E-|-As with unreliable bulge-disk decompositions. 
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TABLE 6 
Color Gradients 





d(B-R) 
dlog r 




d(B-R) 1 
dlogr ]■ 

(mag dex~ ) 


lout 

(mag dex ) 




ID 


/ n 1 \ 

(mag dex ) 


(kpc) 


Morphology* 


EAOIB 


13 + ni 








p 


EA02 


no 4- {Y\ 

U.UU IE. u.uo 


u . uo _I_ u . uo 


— n 9^ + 04 


"^i^ -1- 


N,D 


EA03 


n on n 01 








p 


EA04 


_n on -k n m 


fi ni -1- 

u. u ± ^ u . uo 


n n? + 01 

u.u ( in u. u -L 


_n 90 -1- 0'^ 


N D 


EA05 


ZxL U.UO 


U. O / _U U . 


n Qf^ + n 07 
u. yo ^ u. u / 


n 97 -1- n 09 
u. ^ / u. u^ 


P BC 


EA06 


n 44 + n f)'^ 


n 1 Q _|_ A no 

U.-LiJ _1_ L/.Uii 


U.UU _1_ U.UiJ 


_n 05 ± 05 

U.UU _L U.UU 


P BC D 


EA07 




n 9« -u 0*^ 


—0 4fi -1- 03 

— U.^U ZL U.UO 


_n 09 -1- 01 

— U.U^ ^ U.UX 


N D 


EA08 


n m -1- n m 


U.U^ ^ U.UO 


_n 04 -1- 01 

— U.U1 ZL U.UX 


n 40 _|_ n 1 1 


p 


EA09 


n 06 + n 04 


—0 .38 + 01 


1 02 + f) 05 


—0 18 + 01 


P,BC 


EAIO 


0.08 ±0.01 


0.66 ± 0.03 


0.06 ±0.01 


0.63 ± 0.09 


'p 


EAll 


0.17 ±0.04 


0.13 ± 0.03 


-0.09 ± 0.03 


0.43 ± 0.03 


F 


EA12 


-0.31 ±0.05 


0.08 ± 0.02 


0.34 ± 0.04 


-0.68 ± 0.03 


N,D 


EA13 


0.08 ±0.02 


-0.11 ±0.03 


0.33 ± 0.02 


-0.02 ± 0.01 


F,D 
P,BC 


EA14 


0.29 ±0.03 


0.13 ±0.01 


0.75 ± 0.02 


0.08 ± 0.02 


EA15 


0.12 ±0.01 








P 


EA16 


0.45 ±0.02 


0.07 ±0.02 


0.63 ±0.02 


0.28 ±0.02 


P,BC 


EA17 


0.18 ±0.02 


-0.43 ± 0.03 


0.72 ±0.06 


0.10 ±0.01 


P,BC,D 


EA18 


0.06 ±0.02 


-0.15 ±0.04 


0.40 ±0.04 


-0.10 ±0.02 


F,D 


EA19 


-0.25 ± 0.02 


0.46 ± 0.03 


-0.37 ±0.02 


0.23 ± 0.06 


N,D 


EA21 


0.09 ±0.01 


0.13 ±0.09 


0.15 ±0.02 


0.04 ±0.02 


F 


Note. 


— Because 


a single linear fit is not the best representation of most 



color profiles, the statistical uncertainties are always too small and unrealistic. 
Therefore, we calculate bootstrap uncertainties. 

Morphology of the color profile: N (negative), P (positive), BC (blue-core), 
F(flat or variable), D(dust) 
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APPENDIX 

DETAILED MORPHOLOGIES OF E+A GALAXIES 

High resolution HST images enable us to identify a wealth of small and large scale features. For example, the low 
surface brightness tidal features, small companion galaxies, blue cores, bars, and even point-like, newly-formed star 
clusters are essential clues in understanding the causes and end-products of the E+A phase. In this appendix, we 
describe the qualitative morphologies of 21 E-t-A galaxies. Given that our E-t-A sample was selected using uniform 
spectroscopic criteria (large Balmer absorption and no [O II] emission), it is striking that the morphologies of E+A 
galaxies are so diverse. 

Tidal Features and Disturbed Morphologies 

Eight E+As (EAOlAB, 02, 03, 07, 08, 10, 12, 13) show clear tidal features indicative of recent galaxy-galaxy 
interactions. These features include tidal tails that extend over a few tens of kpc and bridges that connect apparent 
companions. In addition to these dramatic large-scale tidal features, another three (EA04, 11, 18) exhibit disturbed 
morphologies such as shell-like structures (EA04 and EAll) and highly irregular isophotes (EA18). We conclude that 
the fraction of E+A galaxies with readily identifiable, i.e., brighter than our detection limit (fiR < 25.1 ± 0.5 mag 
arcsec" ^), merger/interaction s ignatures is ~ 55% (±15). 

Since IZabludoff et al.l (|1996f) first claimed galaxy-galaxy interactions/mergers as the main mechanis m for E+A 
forma tion, other studies using larger samples of E+As have reached similar conclusions. For example, iBlake et al.l 
(|2004') found that -13% of E+As selected fr om the 2dFGRS, us ing th e same spectroscopic criteria as was done for 
the LCRS sample, show tidal features or disturbed morphologies. iGotol (|2005D found that ^--^30% of his sample shows 
tidal features. The higher merger fraction that we report here is almost certainly due to the improved sensitivity and 
higher resolution of HST imaging rather than intrinsic differences among E+A samples. For example, the low surface 
brightness features seen in EA08 and the shell-like structure in EA04 would not be detected in typical ground-based 
imaging. 

Do merger features correlate with increasing EA numbers (i.e., increasing D^QQa), which is likely to correlate with 
the mean stellar age of the galaxy? We find tidal features more often among E+As with smaller EA numbers. The 
distribution of EA numbers for E+As with tidal features is not random at the ~90% confidence level. Because the 
duration of the E+A phase (~ Gyr) is longer than the dynamical relaxation time of the merger, we expect this 
correlation if EA numbers constitute a rough age sequence. 

Companion Galaxies 

Among the 11 E+As with interaction/merger signatures, five (EAOlAB, 10, 11, 12, 13) have companion galaxies 
within ~ 30 kpc that appear to be interacting with the E+A. The EAOlAB system is spectroscopically confirmed 
(|Yang et al.l 1200(1 ) . but the others may only be projected companions. For example, EAl5's companion (R.A. = 
14:40:45.2, DEC. = —06:39:53), which we also observed in our slit, is at a different redshift. The properties of the 
projected companion galaxies are as diverse as those of the E+As. They range from almost 1:3 mergers (EAOlAB) 
to apparent 1:75 minor mergers (EAIO) as judged from the relative R band fluxes. Are the interactions with these 
companion galaxies responsible for the starburst in E+As that ended < Gyr ago? If these apparent on-going interac- 
tions produce the E+A spectra, then there was a close interaction < Gyr ago, and we are now witnessing a subsequent 
passage. Therefore, a wide range of merger configurations may cause the E+A phenomenon. Alternatively, another 
interaction/trigger produces the E+A spectra, and the current interactions are solely by chance. Larger surveys in- 
cluding redshift measurements of the companion galaxies will help to resolve this issue. Conversely, not all interactions 
lead to E+A's. Yagi et al. (2006) present a spectroscopically confirmed interacting E+A system where one galaxy is 
an E+A, but the other has neither current star formation nor post-starburst signatures. 

Dust Features 

Seven E+As (EA06, 07, 12, 13, 17, 18, 19) show dust features, such as lanes and filamentary structures, in the 
two-color composites or the residual images. In the WFPC2 sample, we also suspect that EA02 and EA04 might have 
dust from their irregular residual images, but the shallow B439 images prohibit us from confirming it. Signatures of 
dust, although typically fairly minor, are present in over a third (33-43%) of our E+As. 

Th ere have been sug gestions that star-forming galaxies could be disguised as E+A galaxies due to obscuration by 
dust (ISmail et al.lll993 ). However, this is not an issue for the E+As in this sample. As discussed in 21 only three E+As 
(EA07, 12, 19) have color profiles seriously affected by dust lanes [A{B — i?) > 0.5]. Furthermore, radio continuum 
emission is detected from only t wo of the dusty E+As (EA12 and 19) with inferred star formation rates (SFR) of 5.9 
and 2.2 Mq yr"-'^, respectively ()Miller fc OwenI [20011 ). EA07 was not observed in the radio. We conclude that the 
contamination rate due to the E+A selection criteria ado pted b y Zabludoff et al. (1996), which employs the average 
Balmer line strength and [O II], is less than 15% (see also iBlake et al.ll2004l : IGotoll200'l . 

Barred Galaxies 

Only three E+As galaxies (EA05, 14, 21) exhibit clear signs of bars in the residual images, while two other E+As 
(EAll and EA19) have elongated residuals that suggest the presence of a bar. Because there is no model profile for 
the bar component, we fit these galaxies as well as possible using the various model GALFIT components. Because 
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introducing too many free components makes it impossible to interpret the components physically (|Peng et al.ll200^ . 
we use only three-component models, each with a combination of two Sersic profiles and one exponential disk profile. 
We show the best-fit model profiles in Figure [3] for the barred E-l-A galaxies. The bars appear to be fit reasonably 
well with a steep (n < 0.8) and flat (g ~ 0.5) component. 

POST-BURST AGE DETERMINATION 

The ages of the young star cluster populations in E+A galaxies could be used to break the degeneracy between 
the burst strength and the time elapsed since the starburst, because a simple stellar population (SSP) is an excellent 
approximation for the star clusters and therefore no additional assumptions are required to model the underlying old 
stellar populations. To derive the time since the starburst (post-burst age), when the majority of the detected clusters 
and young stars formed, we use a statistical test and the simple assumption that all star clusters formed during a 
single burst. Therefore, we assume that their spread in colors is due entirely to color measurement errors. For each 
post-burst age, we draw a thousand sets of N cluster colors, each scattered according to a Gaussian error distribution 
described by the measurement error cr*^^, where N is the number of cluster candidates found in each E-f A. These 
simulated colors are compared to the distributions of observed colors. We determine the fraction of sets that are 
statistically indistinguishable from the observations using the two sample K-S test. If the fraction of acceptable sets 
is less than 5%, we reject that model age as the time since the starburst. 

Next we test whether the key assumptions in determining the post-burst ages — 1) the Gaussian distribution of 
observed colors and 2) the single age cluster population — are statistically acceptable. First, we consider whether 
the assumption that the measured cluster colors are randomly distributed according to a Gaussian is realistic. Due 
to small number statistics in each galaxy, it is not possible to test this assumption for the clusters in each galaxy. 
Therefore, we apply the following test to all clusters discovered in the five E-|-As (EA07, 08, 11, 15, 18). For each 
E-l-A, we calculate the deviation of each cluster color from the median color of the sample in terms of its measurement 
error. Si — {Ci — Cmed)/<^lrn where Ci and Cmed are the measured color of the cluster and the median color of clusters 
in the given E-|-A, and dg^^ is the error in each cluster color including uncertainties in background subtraction. We 
examine the distribution of Si's. If the Si's follow a Gaussian distribution with zero mean and a standard deviation 
of unity, our assumption of Gaussian distributed errors can be justified. If not, e.g., if the Si's show a skewed or 
flat distribution, their measured colors could be seriously affected by the other factors. In Figure [191 we show the 
distribution of Si's. The histogram is marginally consistent with a Gaussian distribution with unit standard deviation 
(cTfit = 1-34), therefore we conclude that there is no strong evidence against this assumption. 

Second, we test the validity of the assumption that all of the detected star clusters formed during a single instan- 
taneous burst. We generate a thousand sets of simulated cluster colors with an age spread of At. Gaussian random 
errors with a dispersion corresponding to the measurement errors al„ are added to the simulated colors. Now we 
calculate the Si 's using these simulated sets of cluster colors and compare the width afn of the Si distribution with the 
observed value am = 1.34. If the observed aat is signiflcantly different from the distribution of CTfit, then we reject the 
given At. We flnd that Ai < 0.1 Gyr is rejected at the 95% confidence level. However, the small number statistics and 
low S/N of the color measurements do not allow us to constrain the upper limits of the burst durations. Therefore, 
the instantaneous burst assumption might not be strictly true, and the burst duration should be at least longer than 
100 Myr. On the other hand, using an E-l-A galaxy sample selected from the Sloan Digital Sky Survey, lYang et"al] 
([2008i) show that the timescale over which star formation ends tends to be smaller than ~ 200 Myr, and larger burst 
durations (several hundred Myr) would wipe out the correlations that we find among post-burst ages, global E-l-A 
colors (i.e., EA number), and the bright end of the star cluster luminosity function. Therefore, we argue that the 
single burst age approximation is valid to within a few hundred Myr. 
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Fig. 19. — The composite distribution of Si = {d — Cmed) I o'lrr from five ACS E+As with cluster candidates. The histogram is marginally 
consistent with a Gaussian with cr = 1. There is no strong indication against our assumption of Gaussian-distributed uncertainties. 



